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SOME INTEGRAL INEQUALITIES CONNECTED / 
WITH THE CALCULUS OF VARIATIONS 


By G. H. HARDY AND J. E. LITTLEWOOD (Cambridge) 
[Received 15 July 1932] 


1. THE inequalities which we consider in this note may appear a 
rather miscellaneous collection; but they have this in common, that, 
unlike most of the inequalities of ordinary analysis, they depend 
upon, or are suggested by, problems in the calculus of variations. 
In this they resemble the inequality discussed by Bliss and by our- 
selves in vol. v of the Journal of the London Mathematical Society;+ 
but the variational theory underlying this last inequality is decidedly 
more abstruse than any to which we shall appeal. The simpler 
problems discussed here should help to elucidate the difficult argu- 
ment of Bliss. 

Our proofs will be of three types, which we call A, B, C. Proofs 
of type A are strictly elementary; however the results may have 
been suggested, the proofs do not depend on variational theory. 
Proofs of type B depend on our knowledge of the form of an extremal, 
i.e. of a solution of the Euler equation associated with a variational 
problem, but are elementary in all other respects. It is natural, in 
view of the ‘genesis of Euler’s equation, that the principal weapon 
used in such proofs should be partial integration. Finally, proofs of 
class C depend essentially on the ideas and methods proper to the 
calculus of variations (as does the proof of Bliss to which we have 
referred). The interest of the paper lies largely in the logical relations 
between these different types of proof. 

In order to avoid tiresome repetitions we state, once for all, that 
whenever y and its derivative y'(or f and f’) occur in a theorem or its 
proof, it is assumed that y is the integral of y’ (or f of f’). Similarly, 
if y” oceurs, y’ is the integral of y”. This assumption is essential in 
all such variational problems. 

We shall require the following lemma: if y’ is L? in (0, a), then 
y = 0(x!) for small x; and if a = 00, then y = o(x!) also for large x. 


+ Bliss (3), Hardy and Littlewood (6). 
t y’ is measurable and its square integrable. 
R 
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In the first place, by the inequality of Schwarz, 


=(fywal. < «| ytd = 00 v) 
“g 
for small x. Secondly, if a = oo, we can choose X = X(e) so that 
{ y’? dt - 
x 
and then 


fy(x)—y(X)}? = Ja [ y’ tat) < (x—X) | y dt < ex, 
xX 


ly(a)| < |y(X)|+ea! < eat, 
for x > a(e,X) = 2 (e). 
2. We begin by proving 
THEOREM 1. If w > 4, y(0) = 0, y(1) = 1, and y’ is L?, then 
1 


J = Jy) = [ (m4 


where 


There is equality if and only if y 
Our proof is of type B. We begin - supposing p > 4, so that 
a> 0. The Euler equation for J is 


parry" +-y = 0, 
and y= = at 
is the extremal satisfying the conditions. There is equality in (2.1) 
when y = 

We cies ge erally y= 
so that 2(0) = 2(1) = 0. Then 
J(y) = J(n)-+F(2)-+2K(n,2), 
1 


where K(n,z) = [ {un’s’— 23) dx. 
a 
0 


Since 2’ is L?, z = o(x') for small x. Hence, integrating by parts and 
observing that 7 satisfies (2.2), we find that K = 0 and 


J(y) = J(n) +I). (2.3) 
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It is therefore sufficient to prove that 
J(z) >0 (2.4) 
unless z = 
We now write 2 =. 9u, 
so that u(1) = 0. Then 


mt 


od 


1 1 1 
J;(z) = | (ue) dx = p | (nu’+-n’u)? dx — | dx 
C) C) 


C) 
1 1 


1 
2 
=~ | Pu’? dx + | (un’*—) u? dx + 2p | nun’ dx. 
3 5 





But 1 1 
2 | nnn’ dx = —p(n7'U*)s—p | (92+ 9" )u® de. 
8 r) 
Combining these equations, and observing that 
‘= (ha), "+t =0 uw’ = —1t4, 
n = (3+ i. eyt[a=o 9 ae 
and z = o(2'), we obtain . 
2 rye 
Ie) = —nd+ay(Z) +o | (m’) ae 
3 


1 1 
J(z) = | (ue—5) up | (v— 2%.) ae. (2.5) 
0 0 


This is positive unless the integrand is nul, when z = Cn, which is 
impossible unless C = 0. This proves (2.4), and therefore the 
theorem, when p > 4. 

It will be observed that (2.4) and (2.5) give 


1 1 

9 ‘ 2 , 3 a 5 

J(y) = | (uy —) dx = ate | (v —2* v) dx. (2.6) 
0 0 


This identity may naturally be verified directly by partial integra- 
tion when its form is observed. Since 7 has now disappeared, we 
may take pp = 4, a = 0 in (2.5) and (2.6); and this completes the 
proof. 

Trivial transformations of (2.6) give 


é é 
A te oe oe ,_ ite)? 2.7 
0 


0 
R2 
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where now y(0) = 0, y(€) = c. Suppose now that y’ is L? in (0,00), 
so that c = y(é) = o(€) for large €. Then, if we make € > 0, (2.7) 


becomes = a 
= y? A ee ; 
y’2— 2) dx = of ee dx. 
F (otha mn fp tess) 
0 0 


Taking » = 4, we obtain 


ft ae y 
dy?) dx - ‘| — 4) dz, 
| ( d “) J (v zy 
0 0 


which is positive unless y is nul. This identity} gives the simplest 


proof of the inequality 


dx < 4 y"? dx. (2.9) 
a 
0 0 

3. Our proof of Theorem 1 is of type B; underlying it, there is 
a proof of type C. The family of extremals 

y= ay xtra 

defines a field{ surrounding the extremal y = 7. The slope-function 
of the field is 
j yy 


9 = a($-+a)x-t+4 — (4-+a)= 
‘s r 2 


and the Weierstrass excess-function is 
E = p(y’—p)?. 


Hilbert’s invariant integral is 


J*(y) = ial {(F —pF,) dx +F, dy}, 


where 
and the general theory tells us that 

J(y)—I(n) = J(y)—J*(n) = J(y)—JI*y) = [ € de, 
which is (2.6). It will be found that in this case 


(F—pF,) dx +F, dy = a E ) 
(3 —a)a 


{ For similar, but less simple, identities for series, see Grandjot (5). 
t See the treatises of Bolza (4) and Bliss (2). The field is ‘uneigentlich’, 
since all the curves pass through the origin. 
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4. Similar analysis may be applied to the more general inequality 


jefe <(ehy frm 
0 0 


where k > 1, y’ > 0 


THEOREM 2. If 


ee | k \* & 
* nig a 1 > —_ 
k>1, ETP = Il, hw > (=) os XZ, 


y(0) = 0, y(1) = 1, and y’ is positive and L*, then 


J ag dx > . 
w= | (1-5) @ > gasp 
0 
where X is the (unique) number defined by 
p(k—1)A*-1(A— 1) +1 = 0, A> lk’. 
The underlying identity, corresponding to (2.6), is 


1 yk . i 
| (w*-S) = gaat 


0 


on ora Ye 


and the integral on the right-hand side is positive, unless y’ = Ay/x. 
When p = K, A= Ik’ and pA¥ = 1. 


5. Tororem 3. If y(0) = y(1) = 0 and y’ is L*, then 


1 
faatgee st | ve 
0 


with equality only when y = Cx(1—2). 


: 2 
Here J(y) = (wv2- = z ) dx 
J x(1—z) 


and Euler’s equation is 9 
y+—*_ = 
a(1—z) 


The general solution is 


y= x(1 —2)(0-+ Blog a) 
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and the family y = ax(1—z) 

defines a field} of extremals round any one of its members. 
Here — 

a(1—z) y 


p = a(1—2z) 


" —_ = 
(F—pF,) dx +F, dy = al —- r’) 


2a( 1—2z) 


and the identity to which we are led is 


1 
1. ee ) dx = : (v PR: ) dx. 
a(1—z) 2. a(1—z2) 


Our proof is of type C or B, according as we construct this identity 
from the general theory or verify it directly. 
Theorem 3 was suggested to us by the following theorem of Pélya.t 


9 


THEOREM 4. Jf 2(1—«x)y’? is integrable, then 
ES 1 2 : 
0 < | y? dx -(| y de) <} | a(1—ax)y’? dx. 
0 0 0 
The sign of equality can occur in the two possible places only when y is 
(i) constant, (ii) linear. 
Pélya’s proof is as follows. The first inequality follows at once 
from Schwarz’s inequality. For the second, we have 
1 1 2 1} 
| y* dx —(| y dz] = 4% | | (y(u)—y(v))* dudv 
0 0 0 
> -4 2 . 3 (3 
= | du | dv (| ya < | du (v—u) dv y’? dt 
0 u u 0 u u 
4 3 ; 
= | y"? dt | du | (v—u) dv = } t(1—t)y’? dt, 
0 0 t 0 
with equality only when y’ is constant. 

Theorems 3 and 4 can also be deduced from ‘Parseval’s theorem’ 
in the theory of Legendre series. A trivial transformation reduces 
them to the forms , 

| es dx < } | y"* dx, (5.1) 
-1 -1 


+ Here ‘uneigentlich’ at both ends. t Communicated to us by letter. 











SOME INTEGRAL INEQUALITIES 
with y(—1) = y(1) = 0; and 


1 1 
[ yrdx <4 | (1—2*)y’? dx, 
~4 ~—% 
j 
with ydx = 0. 
—% 
To prove (5.1), we suppose that 
y ~ >a, P(x 
with ag = 0. Then 
x 
y= >a, | PW dt = > a,p,(2), 
5 
say. It can easily be verified that 
2 
PmPn lx = 0 : _ 4 i te 
| 1—a sabi. n(n-+-1)(2n-+1) 
—1 


and (5.1) then follows from 


wo 


c 2 az a 
Pe: a n 
3 n(n+1) 2n+1 ~ hi 2 


1 
Similarly, if we suppose that 
- be 2 > ay, P,, (x 
and use the formulae 
, _ 2n(n+1) 
[ I—2*)P_, Pade =0 (m An), = 
= 


(5.2) will be found to follow from 
ro a < > n(n-+-1) _& 
2n+1 2 2n+1° 
In each case there is equality only when a, = a, = ... = 0. 


6. TuroreM 5. Jf y(0) = 0, y(1) = 1, k is a positive integer, and 
y’ is L?*, then 


fy dx < o fy dx, 


0 


} ——_ o 
where ht =) ° 


There is equality only for a certain hyperelliptic y. 











248 G. H. HARDY AND J. E. LITTLEWOOD 

If a = 1/C and R 

J(y) = ly ‘2k —ay**) da, 
) 


( 
Euler’s equation is 
(2k—1)y’**-2y” + ay**-1 — 0, 
which gives (2k—1)y’* = b—ay”*, 


If we take b = a, we obtain the special family of extremals 


2h —1\ 12k dy 
“= : 
a) J (yy 


1 1 
dy ] lox. tak 
, a 1—y)-N2kyl2k-1 dy —  cosec—_. 
ea =z | | u) u du 55 OO ar: 
0 


and then, since 


the extremal through the origin will also pass through (1, 1). 
If we denote this extremal by y = », and observe that y’(1) = 0, 
we obtain 


1 


J(n) = [n7’ ak—1) — { n{(2k—1)y ‘2k-2y)" + an*k-l} dx = 0, 


0 
and it is plain that the theorem will follow if y = 7» gives a strong 
minimum for J(y). That this is so follows from the general theory. 
The graph of 7 is a rising curve whose shape resembles generally 
that of the curve y = sin $72, to which it reduces when k = 1. The 
family 
y = an 
defines a field, and the Weierstrass function 
79 > 
fF — hy —p)p* 

is positive; and this is what is required for the application of the 
general theory. 

This proof is of type C and (in view of the trouble of calculating 
the slope-function p in terms of x and y) it might be difficult to 
construct a much more elementary proof. When k = 1, the in- 


equality reduces to 
1 1 


| wae < as [ v2 dx, 
7 . 


0 0 
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with equality for y = sin}. This inequality is very like one due 
to Wirtinger, which has an important application to the proof of the 
isoperimetrical property of the circle.t In this case there are more 
elementary proofs, and the reader will find it instructive to construct 
proofs similar to those of the previous sections, and to deduce the 
result from the theory of Fourier series. 


7. In § 6 we used a proof of type C which it might be difficult to 
reduce to type B. In this section we use a proof of type B which 
we cannot replace by one of type A. 


THEOREM 6. If y and y” are L*, then 
io a) 2 i 3) wo 

(fy da) < a[y dx | y”? dx, (7.1) 
0 0 0 


0 
with equality only when 
y = n = Ce-* SY sin(x sin y—y), (y = 4p). (7.2) 
The Euler equation for the ‘isoperimetrical’ problem of making 


J, = | y’? dx 


a maximum, when 


J,= fy’ dx 


0 


are fixed, is of the form 
ay” +by” +cy — 0, 
and its solutions are linear combinations of real or complex ex- 


ponentials. When we tried to choose the parameters of the solutions 
iu the most advantageous way, we were led to select the function 


(7.2). 


We shall deduce Theorem 6 from 


TuroreM 7. If the conditions of Theorem 6 are satisfied, then y’ is 
L? and a 


Jy) = f y2+y"—y") dx > 0, 


0 


with equality only for y = 7. 


+ See Blaschke (1), 105. 
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In the first place 
x = 
y= dx = [yy — yy” dx. (7.3) 
0 0 
Since J, and J, are finite, the last integral tends to a finite limit when 
X +o. Hence, if J, were infinite, yy’, and a fortiori 4y? = f yy’ dx, 
would tend to +00 with x, which is inconsistent with the con- 
vergence of J). Thus J, is finite, and all three terms in (7.3) tend 
to limits. In particular yy’ tends to a limit, which can only be 
zero. 
It is easy to verify that 
n+7'+n” = 0, = —7"(0) = 0, 
and 


We now write 


choosing ¢ so that 2(0) = 0. Since z differs from y by a function 
which, with its derivatives, vanishes exponentially at infinity, 

zz’ > 0, (n+ 7’ )z > 0, (n+7')z’ > 0. 
Also J (y) = ce? J(n)+2cK(n, z)+J(z), 


and 2 
x: | (nz-+ "2" —n’z') dx 
_ [ (n’+")z dx — [ (n-+-7’)z” dx — [ nz’ dx 
0 0 0 
(n+7'+7")z’ dx = 0. 
Hence (7.4) reduces to J(y) = J(z), and it is sufficient to prove that 
J(z) > 0 


unless z = 0. But, since z(0) = 0 and zz’ > 0, 
i 


. zz" dx, 


and so 2 22+-22"+2"2) dx > 0. 
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This proves Theorem 7. To deduce Theorem 6, we replace x by 
Ax and y(Ax) by y(x), when we obtain 


MJp—A2J,+J_ = | (Aty2—D2y’2+-y"2) dx > 0. 
0 


Since this is true for all A, J? < 4J)J,. 
Theorem 6 suggests that, if k > 1, 


(J yt az) < Ke) f iyit ae fy de 
0 0 0 

We cannot deal with this inequality completely; it is easy to prove 
K(k) finite for every k, but we cannot determine its best possible 
value, even for k = 4. A theorem of Landau? suggests that the limit 
of K"* when k + o should be 4. 


8. It may be interesting to remark in conclusion that the solution 
of the corresponding problem for the interval (—0oo,00) is simpler 
and quite different. 


THEOREM 8. If y and y” are L* in (—o, 0), then 


([v<( fre [ora 


unless y is nul. The constant | is the best possible. 


For, as in § 7, yy’ > 0 when 2 tends to infinity, and 


is) ao 


| y= dz = — | yy” dx, 


so that the result follows from Schwarz’s inequality. 

To prove the constant the best possible, take first 

y=sinx (|x| < nn), =0 (|z| > nz), 
and then round off the corners at x = --nz so as to make y” con- 
tinuous. We can do this with changes in J), J,, and J, which are 
independent of n. Then all three integrals differ from nz by numbers 
independent of n, and 
J? > (l—e)Jg J, 


for every « and sufficiently large n. 


+ Landau (7). 
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NOTES ON LINEAR EQUATIONS IN INFINITE 
MATRICES 


By P. DIENES (London) 


[Received 16 December 1931] 


1. The problem. Definitions. For given y,, the solution of 


PET: =. (n = 1, 2,..., 00) (1) 


hinges on certain matrix characteristics of the table a,,. The same 
equations are used, in the theory of infinite quadratic forms, for 
the transformation of the variables x, into y, in order to reduce 
yn a 
> a, ;.%, 2, to a diagonal form > 5; y3. 

1=1 


n,k=1 
Again, the generalized A-limit of 2, is that of y,, in (1) whenever 


the latter exists. In this connexion, A = (a,,) being given, the 
problem of the consistency of the various generalizations of limit 
leads to the problem of finding a matrix Q and a diagonal matrix D 


such that A—QDOQ-1. (2) 
In the Heisenberg-Dirac theory of atom mechanics the central 
problem is to find an Q and a D such that 

Q-1A0 = D. (3) 

The determination of Q in (2) and (3) requires, as a rule, the solu- 

tons AX—XD=0. (4) 
Also, the fundamental equation for quantization is 

AX—XA = 1. (5) 


Therefore we thought it desirable to collect the more or less obvious 
results on linear matrix equations, with illustrative examples, as 
an introduction to the theories mentioned above. Some results, 
especially in §§ 4, 5, seem to be new. A list of papers on the subject 
is appended for reference. The point of view adopted, first put for- 
ward in A. C. Dixon} (1902) for certain systems (1), aims at estab- 
lishing facts about the existence or non-existence of solutions and at 
developing methods for explicit formulae without making use of 


+ References are to the bibliography at the end of the paper. 
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infinite determinants. For functions of a matrix variable see D. Hil- 
bert (1909), Riesz, chap. 5, and Dienes (1930). I owe some valuable 
suggestions to Professor A. C. Dixon. 

Derinitions. An infinite matrix is a twofold table 


A = (a,,) (n, & = 1, 2.,...,00) 
of real or complex numbers, with addition and multiplication de- 
fined by 


A+B=(Aytby), ch=(cay), AB=(Saybq), (6) 
i= 


where ¢ is any real or complex number (scalar). The first suffix n 
specifies the row, the second suffix the column. The zero matrix 0 is 
the table a,,;. = 0. ‘The elements 5% of the unit matrix, denoted by 
1 or #, are . 
r , {Vi nFAk 


5 — - 
, llifa—k. (7) 


These two definitions are justified by the obvious equations 

0A = AO= 0, EA = AE= A. 
We also notice that the matrices cH form a sub-algebra which is 
simply isomorphic with complex algebra, so that matrix algebra can 
be considered as an extension of complex algebra. 

The elements a,,, form the (principal) diagonal. If, outside the 
diagonal, a,,,, = 0, then A is a diagonal matrix denoted by D, A, 90, T’. 
If a,; = 0 above (below) the diagonal, then A is a lower (upper) 
semi-matrix. If every row (column) contains a finite number of non- 
zero elements as, for example, in lower (upper) semi-matrices, the 
matrix is row-finite (column-finite). If a,,. = 0 outside a certain finite 
number of rows (columns), then A is row-bounded (column-bounded). 
For example, (a,,,) denotes the matrix whose first column is @,,,, all 
the other elements being 0. 

If AB = 1, then B is a rh. (right-hand) reciprocal of A denoted 
by A-}, and A is a lh. (left-hand) reciprocal of B, denoted by —'B. 
Similarly, if both A and B are different from 0 and if AB = 0, then 
B is a rh. 0-divisor of A, denoted by A®, and A is a lh. 0-divisor of B, 
denoted by °B. 


2. Remarks on matrix algebra. The sum of two matrices 
always exists, and 
A+B= B+A, (A+ B)+C = A+(B+0), 


but products are in general not commutative and AB may not exist 
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even when BA exists. For example, if B = (b,,,), then BA = (b,, a,,) 
exists for arbitrary b,, and A; on the other hand 


(AB) ai = =F ayebin (AB),;, — 0, if k = 1, 


so that AB does not exist when > a,,;6;, is divergent. 
7 


The distributive law 
A(B+C) = AB+AC, (B+C)A = BA+CA 
holds in the sense that, if AB and AC exist, then also A(B+C) 
exists and is equal to AB+AC. But A(B+C) may exist when AB 
and AC do not exist. For example, if 
a, = 1 for every n and k, 
b,, = n?+1, c,, = n-*—1 for every k. 

The product law may not hold even for simple matrices. For 
example, if in D and A, d,,, = 0, 8,,,, = 0 and the other diagonal 
elements are different from 0, then we have DA = (d,4,) = 9, 
though neither D nor A is the zero matrix. 


Multiplication is not necessarily associative. For example, if 
a,;, = 1, ¢,, = 1 for every n and k, and if 


pe (> bis) # > (> bis), 


then (AB)C = ») (x b,j) # A(BC) = > (3 ,))- 
= y=1 
As a numerical example we put 
5 ig ee ee (i—1)2 
, Cae P+j-1 G-1f4+j ° 6—-1P4+5-¥ 
in which case the partial sums are 








a2 
8, = By’ 

so that the two different summations lead to two different sums. 

We shall say that a product is associative if (i) the product exists 
for every succession of the multiplications involved, the order of the 
factors remaining the same, (ii) all the products so obtained are 
equal. If, in a class of matrices, every finite product is associative, 
we say that the class is an associative field F. Diagonal matrices 
form an associative (and commutative) field. 
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Bound of a matrix. If, by a given rule, a number |A| > 0 is 
assigned to every matrix of a given class satisfying the conditions 

(i) |cA| = |e||A|, |#| = 1 (c any complex number), 

(ii) |A+B| < |A|+|B|, 

(iii) |AB| < |A||BI, 
we say that the class is bounded according to the given rule, and 
A| is the corresponding bound or modulus of A. For a similar defini- 
tion see Helly, p. 61. Such a bound is the generalization of absolute 
value (modulus) of complex numbers, which justifies the notation. 
We notice that |a,,| denotes the absolute value of the complex 
number @,,,, so that, in general, |a,,,.| ~ |A]. 

The bound is said to be regular if it also satisfies the conditions 


(iv) IGrk| * : A|, 

(v) |D| =d, 
where d is the finite upper bound of |d,,|. Condition (v) means that, 
if the upper bound d of |{d,,| is finite, every regular bound of the 
corresponding diagonal matrix D is exactly d. 

We notice that 


ie 8) nd 
Zz Cp A,| < 2 Cp ae (1) 
p=1 p=1 , 


whenever the right-hand side exists, and in this case, if A,, Ag.... 
belong to the same associative field F, then > c,, A, also belongs to F. 
Examples. Putting > nsi = M.,, > a,,;.| = N;,, the upper bound 
k=1 n=1 
of M,(N,,), when it is finite, is the K,(K,) bound of A (see Dienes 
(1931), chap. 12). 


2) 

The finite upper bound of | > a,,,%,,2;,|, for sequences x, such 
7 nk“n “Vk n 
n,k=1 


2 . 
that > |z, |? = 1, is the H (Hilbert) bound of A. It follows that 
n=1 
« 2) 
la..,.|2 la..|2 2 
Zz |Qnk| < |A| ? 2 \Ank| < |A| ? 
n c=1 
but the converse is not necessarily true. 
It is readily verified that the K,, K,, and H bounds are regular. 
We finally notice that row-finite, column-finite, K,, K,, and H 
matrices, i.e. the class of matrices having a finite K,, K,, or H bound 
respectively, form so many associative fields. 
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In fact, the product of two row-finite (column-finite) matrices is 
a row-finite (column-finite) matrix, and, in 
is) oo io a] 2) ‘ 
> (> a; Dis)Cjx = 2 a,( > bi; Ci)» (2) 
g=1 i=1 / i=1 j=1 
the infinite sums reduce to finite sums. 
If A, B, C are K, matrices, say, 


Oni U 


Ms 


Pin| <|BIY | dy! < [AI |B 


1 


v 
and 


nx 


2 | p> A; 55; ) ex < |C| 2 pH ay ¢b4;| < |C|| B| 2, |@ni! < |A||Bi|C}, 
j i j=11= i= 
so that (2) is absolutely convergent. 

For H matrices see Hilbert (1912), p. 129. 


3. Reciprocals of matrices. The second of the equations 
AX = B, XA=B (1) 


may be written in the form A*X* = B*, where the star indicates 
the transposed matrix (rows and columns exchanged); we therefore 
deal with the first equation only. 

If a diagonal element in D is zero, DX = 1, ie. d,,x,;, = 5% has 
no solution, since, if d,, = 0, say, d,,x,,, = 1 has none. Thus (1) 
may have no solution. However, if A-! is a solution of 


AX = 1, (2) 


A-'B is that of (1) provided AA-B is associative. Therefore we 
shall first examine (2). 

If A-! is a solution of (2), A-!+ A® is another solution, and every 
solution is of that form. Moreover, in the same associative field, no 
A has an A- and an °A, or an —14 and an A®, since (°AA)A-! = 0 
and °4(AA-1) = °A. 

The relation between lh. and rh. ‘reciprocals is illustrated by 
I. (i) If Ais unique and AA-'.A is associative, then A~ is also a th. 
reciprocal of A and the only one for which -1AAA~ is associative. (ii) If 
A has both an A-! and an -1A, and -1AAA~ is associative, then 

14 — A-1, and A has no other two-sided reciprocal for which -\AAA- 
is associative. 

Proor. (i) IfQ isa solution of AX = B, all the other solutions are 
of the form Q+ A®. In particular, the solutions of AX = A are all of 

3695.3 Ss 
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the form 1+ A®. If A-!is unique, there is no A® (since A~1+ A® would 
be another rh. reciprocal), and thus X = 1 is the only solution of 
AX = A. Now, if AA—.A is associative, (A A-1)A = A(A-1.A) = A, 
so that A~!. A = 1. Moreover, if Bisa lh. reciprocal such that BAA 
is associative, B = B(AA-) = (BA)A1 = A-, 

(ii) From -14A'= 1 and AA~! 1, we get 

(7A44)4% = A-, ‘A(AA) = -14, 

so that, when -14AA-! is associative, -14 = A-!. Fixing A-, say, 
every —!A is equal to this fixed A~', i.e. there is one —1A only, and 
similarly for A~!. It may happen, though, that A has no A-! and 
has infinitely many —1A. 

It follows from I that if, in an associative field F, A has a unique 
A-}, this A-1 is a two-sided reciprocal of A, and unique in F. Or, if, 
in an F’, A has both an A~ and an ~14, they coincide and form the 
unique (two-sided) reciprocal of A in F. 


If A has an A- (-1A), and if A is a matrix obtained by ex- 
changing certain rows (columns) in A, the ea exchange of 
columns in A- (rows in -1A) leads to an A- (1A). If A has an A® 
(°A), the same matrix is an A® (°A). 


6; has a solution. Exchange 


ik 


Proor. By hypothesis, ¥ a,,;a 

i 
the rows of A by renumbering them, i.e. by establishing a one-one 
correspondence between n and v, say, indicated by v = f(n), n = g(v) 
from » to v and from v to » respectively. The renumbered equa- 


(3) 


tious are . 
ah 
> Aj(nyi % ik Oj, as (Anni) A ° 

5 i 


Now, renumber the unknowns in 
AX’ = 1, i.e. da 


i 


Mndi Vik ko 
by rewriting these equations in the form 
Ha Ayn)i Vighk) = GK) (4) 
t 


When f(n) = &, then n = g(k) and thus, for 2/4, = x,;,, (4) coincides 
with (3), i.e. xq) = x, is a solution of AX’ = 1. 

This result extends to AX = B in the form: if (1) has a solution 
viz, then 2,4.) is a solution of (a,,),) X = (bypyu9)- 

The second statement in II follows from the fact that, in AX = 0, 
every column 2;;,, 7 = 1, 2,...,00, satisfies the same equations. 

We also notice that if A is a matrix obtained by ruling out certain 
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rows of A, every solution of (1) is a solution of AX = B, where B is 
B without the rows corresponding to the missing rows in A. In 
particular, if A has an A-1, A has an A-, and if A has an -1A, then, 
leaving out columns in A, we obtain matrices having a lh. reciprocal. 
It follows, for instance, that if A has no A-', no addition of new 
rows leads to a matrix with a rh. reciprocal. 


Ill. If A and B have rh. (lh.) reciprocals such that ABB-1.A-1 
(1B.-1AAB) is associative, B-!. A (-1B.-1A) is arh. (lh.) reciprocal - 
of AB. If there is an 1A and a (AB)- such that AA B(AB)-.A 
is associative, (A B)-!.A isa B-. 

Proor. (AB)B-!.A-1= A(BB)A+=1. If now we multiply 
A B(AB)-! = 1 by -1A on the left and by A on the right, we obtain 
B[(AB)-, A)] = 1. 

We notice that when A and B belong to the same associative field 
F, and A-! and B-' are unique on the right, B-!. A- is the unique 
(A B)- in F. In fact, in this case, A~! and B- are also lh. reciprocals 
and thus it follows from (A B)X = 1 that X = 7B.-14. 


IV. If A and (14+ BA-)- exist and (A+ B)A-\(14+ BA-)* is 
associative, A-1.(1+-BA-)-1 is an (A+B). If -1A, -(14+-1AB) 
exist, and —\(1+-1A B).-1A(A-+ B) is associative, -1(1+--1A B).-1A is 
an (A+B). 

In these statements A and B may be exchanged. 

Proor. (4+ B)A-1.(14+ BA)! = (14 BA4)(14+- BA)1 = 1. 

We notice that if A-1, B-, (1+AB-)-1, and (1+ BA-")- are 
two-sided in an F, then 

Bo. (14AB>)> = A. (14+ BA4)4, 
In fact, from (1+AB-)B = (14+ BA-)A = A+B, we have 
1+ABoO= (1+ BA)AB, (14+ BA-)-1.(1+AB") = AB, 
A-1,(1+ BA-!)-1,(1+AB-!) = B-, ete. 


We notice that exchanging the first two rows is the same as multi- 
plying A on the left by the matrix 


0 1 0 0 
1 0 0 0 
0 0 1 0 
0 0 0 1 


$2 
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and we rule out the second row, say, of A if we multiply it by the 


matrix 2 60 0 2 
et oe Ge | eres 
6 8 0 1 


Ly 
In this way II can be established, at least in simple cases, by III. 
I owe this remark to Professor A. C. Dixon. 


4. Examples. Applications. We shall now give some matrices 
having no reciprocal, and some having some kind of reciprocal. 

If A has a 0 row (column), A has no rh. (lh.) reciprocal, since 
> 4, %, = df reduces to 0 = 1 when k =n’ for which a,,; = 0, 
$= 32 00. 

If the terms of a row (column) in A are the same multiple of another 
row (column), A has no rh. (lh.) reciprocal, since if ag; = ca,),, say, 
> a9; 2, = 5% becomes cdi, = 87, i.e. 0 = 1 when k = 2. 

, 


D 
More generally, if } c;a;; = 0 and c, ¥ 0, say, then > a), 2; = 3; 
71 i 


ir” ; ; 
becomes — — S c; 5}, = 8}, ie. 0 = 1, when k = 1, provided the two 
Cy — " 
J =2 


wn eo) 
summations in } > c;4;;,%,,, can be interchanged. This is the case 
i=1 9=2 
if the second sum is finite, or if c; considered as a diagonal matrix, a;, 
and 2:,; all belong to the same associative field. In an associative field 
with a regular bound, this condition only means that |c;| is bounded. 

We say that A is a cross-matrix, if all its non-vanishing terms are 
in a finite number of rows and columns, as, for instance, row- 
bounded or column-bounded matrices. 

No cross-matrix has a reciprocal, 

Proor. By hypothesis, there is an r such that @,,,,.,;= 0 for 
every n and k. Take r+1 rows such that a,; = 0 for 71> r, and 
take them for the first r+] rows. Let A denote the determinant (q@,,;.) 
(n,k = 1,2,...,r). If A = 0, there is a linear relation between the 

r , 
r linear forms > a,j, (n = 1, 2,...,7), leading to > c, 5). = 0, which 
i=1 i=1 


is contradictory for k = j if c; ~ 0 (and there is such a ¢; unless all 
= 

the a, vanish). If A + 0, the equations ¥ a, ; x, = 3 (n = 1, 2.,...,7) 
i=1 


are equivalent to A.z,; = A,,, where A,, is A in which the nth 
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column has been replaced by 8% (Cramer’s rule). Since A,,,,., = 0, 


r 
it follows that 2,,,;=0 for i>0, and thus }a,,,;%), = 67 
m1 


reduces to 0 = 1 fork = r+1. 
In case of row-finite matrices, see Toeplitz (1909), AX = 1 re- 


duces to ‘a 
da -2,, = 8h 
a  niXix he? 
=1 


so that in many cases we can solve these equations one by one or 
in groups. For example, if the rows are of unequal length and there 
is no zero row, we can arrange them in increasing order to form 
a stair-like matrix. We can also dispose of the arbitrary terms to 
obtain a row-finite reciprocal. 

In particular, for a lower (upper) semi-matrix where a,,,, 4 0, we 
obtain a unique rh. (lh.) reciprocal which is a lower (upper) semi- 
matrix with a,) = l/a It follows that this rh. (lh.) reciprocal is 
at the same time a lh. (rh.) reciprocal and is the only two-sided 
reciprocal among lower (upper) semi-matrices. 

The same lower (upper) semi-matrix A, however, may still have 
infinitely many lh. (rh.) reciprocals. For example, if the non- 
vanishing terms are @;;_,, @;;, (two diagonals), XA = 1 reduces to 

Ty Ki Ute t lnk Ue = 
which, by fixing n and varying k, determines the nth row of X when 
the first column of X is arbitrarily given. 

If the matrix has no special form leading to simple equations for 
its reciprocal, we have to restrict ourselves to an associative field 


nn* 


wo 
with a bound. We say that > A® is (absolutely) convergent, if 
p=1 
> a) is (absolutely) convergent for every n and k. 
p=1 


I. If A is a matrix of a bounded associative field F and if |A| < 1, 
then 1+ ¥ (—1)?A” is absolutely convergent, its sum is the unique, 


p=1 


and two-sided, reciprocal of 14+-A in F, and 
1 


> ieel << 
2(-» l<i=a 


hs 


For regular bounds such that |A| < 4, we also have 


= 1 
14 $(—1r42|> 7). 
p+ 2 (wre > ra 
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2) 
Proor. It follows from |A”| < |A|” that, if > |c,||A|? = ¢ and 
1 
x e 
> |9,||B\? = g are convergent, C = }c,, A” and G = }g, B? are 
p=l Pp p 
matrices of F such that |C| < c, |G! <g. Moreover, it follows from 


(CG)nrx ‘ ‘ep Cy (A”),5 z Iq (BY), < cq, (3) 
a qd 


that the terms of the product > c, A” > g, B4 can be rearranged at 
will, without changing its sum. Applying this remark to C = 1+A 
x 
and G I+ > (—1)”A?”, we readily prove the first statement. 
9=1 


Also, for regular bounds, it follows from |A| < 1 that 
1+A| > |(1+A),,.| > 1—|a,,; 

Thus, putting B = (1+A)-1, we have 

A| 
1—|A| 
Since B satisfies the equation B = 1—AB, we have 

B| > 1—|AB| > 1—|A]||B| > 0, 

= ] 

1+|A| 


AB| < |A||B) < 


and thus B 


An immediate extension of I, by 3. IV, is 


Il. If, in a bounded associative field F, an A-* exists such that 


|BA-1| < 1, then A-[1+ S$ (- 1)°(BA | is a rh. reciprocal of 


p=l1 


L 


A+Bin F. If A-‘is the only rh. reciprocal of A in F and |BA-| < 1, 
> (—1)(BA)| is the unique, and two-sided, rh. reci- 


then A {1 |. 
p=1 


procal of A+-B in F. 

For H matrices, see Riesz, p. 113. 

If the lower bound of |d,,| is d > 0, D has a unique reciprocal and 
any regular bound of D-! is 1/d. Hence, if |B| < d (regular bound), 
|BD-!| < |B||D-1| < 1 and II proves that D+ B has a unique and 
two-sided reciprocal in F containing both D and B. For matrices 
which are K, and K, matrices at the same time, this result has been 
recently proved by R. G. Cooke (1931). It follows for an associative 
field F with a regular bound that if 0 << d<a,, <d’ and |A®| < d, 
where A® is A in which the diagonal terms have been replaced by 


0, then A has a unique reciprocal in F. 
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Putting x,, = 2, %; = 0 if + > 1, b,, = 6,, 6,, = 0 if + > 1, we 
can write x 
- PT =x 6, (4) 
in the matrix form AX=B8B (5) 
with the proviso that we consider only solutions of (5) in which all 
the terms vanish beyond the first column. It follows that, if A has 
an A~! such that AA~-*. B is associative, A-!. B is a solution of (4), 
since with B also A-!. B is a first-column matrix. Moreover, if A-! 
is unique in an associative field F, A-1. B is the unique solution of 
(4) in F. 
For example, if A is a K, matrix with a K, reciprocal A-1, and if 
the sequence |b,| is bounded, A-!. B exists and 2, = (A-!. B),, is 


“ni 


a solution of (4) such that |2,| is bounded. If A is a K, matrix 
with a K, reciprocal A-', and if > |b,,| is convergent, 2, = (A-!. B),, 
n 


is a solution of (4) such that ¥ |2,| is convergent. 
k 


These remarks apply in particular if A, D, and B belong to an 
associative field F with a regular bound, and if |A —D| is less than the 
lower bound of |d,|. If A is a K, matrix, D = 1, |b,,| bounded, we 
obtain a slight extension of Dixon’s first result (Dixon (1902), p. 191). 

We close this article by some remarks tending to show that, in 
infinite systems of equations symbolized by 

AX = B, (6) 
the transformation of A into a diagonal matrix D = Q-'. AQ corre- 
sponds to the elimination process. 

If A, B, C are in F, then (6) and AXC = BC or CAX = CB are 
equivalent in F, i.e. have the same solutions in F provided that C 
has in F a rh. or a lh. reciprocal respectively. Therefore, if there is 
an Q-' in the field F containing A, B, and Q, in which case Q-! has 
a lh. reciprocal in F, viz. Q, then (6) and 

Q-1,4Q.01.XQ = 2-4. BQ, ie. DY = B (7) 
are equivalent. 

Now, if no d,, vanishes, the only solution of (7) in F is Y = D“.B 
and thus the only solution of (6) in F is 

X = QD-. BO- = QD+.0-1. B. 
If d, vanishes for the suffixes n’, say, the necessary and sufficient 
condition that (7), and thus also (6), have a solution in F is that 
b,,, = 0 for every k, and in this case the rows 2,,,, are arbitrary. 
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Since Q-?, AQ = D is a diagonal matrix, (Q-?.AQ),;, = 0 unless 
k =n. Therefore, if d,, = 0, we have 
(Q-1, AQ),,,, = 0 for every k. : (8) 
On the other hand, the explicit form of the conditions b,,,, = 0 is 
(Q-1. BQ),,,, = 0 for every k, (9) 
i.e. when d,, = 0 for some affixes, A and B have to satisfy the same 
conditions in order that (6) should have a solution. This is the exten- 
sion of a well-known property of finite systems to a class of infinite 
systems of simultaneous equations. 

An obvious generalization of this result is as follows. If (i) A, B, 
Q, W belong to the same associative field F, (ii) Q has a rh., W a bh. 
reciprocal in F, (iii) WAQ = D, then (6) and 

WAQ.Q+, XO = WBA, ie. DY = B, (10) 
are equivalent and the same conclusion holds as before with con- 
ditions (8) and (9) replaced by 

(WAQ),,,=0,  (WBQ),=0, for every k. (11) 


5. Two general types of linear equations. If X is a solution of 
the linear equation AXA’ BXB' =C. (1) 
every solution is of the form X-+-Q, where Q is a solution of 

AXA'+ BX B' = 0. (2) 
Multiplying by B-' on the left and by A’- on the right, (1) and 
(2) assume the simpler forms 
AX+XB= C, (3) 
AX+XB= 0. (4) 
For (4) we have the following results. If Q is a solution of (4) and 
C any matrix commutable with A (AC = CA) such that ACO and 
COB are associative, then CQ is another solution of (4). In fact 
ACQ+ COB = C(AQ+QOB) = 0. 
In a like manner QC is a solution if CB = BC. 

Conversely, if CQ is another solution of (4), ie. if AQN+QB = 0 
and ACQ+COQB = 0 and the products are associative. we have 
ACQ = CAQ. Thus, when there is an Q-!, CA = AC. We finally 
remark that if there is an “1Q, and if W.-1 
be written in the form W = CQ, for W = X 
fore, in that case, all the solutions of (4) are of the form CQ, where 
CA = AC. This is the case in an associative field ¥. Hence 


2 is associative, W can 


¢ 
Q can be solved. There- 
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|. If A and B are in an associative field F and Q is a solution of 

(4) in F possessing a two-sided Q- in F, all the solutions of (4) in F 

are of the form CQ, where C is any matrix in F commutable with A, 
and all such products are solutions of (4). 

Suppose now that A, B, Q are in F. If there is a B-! and an Q-! 
or an Q in F, then —QB-1.Q-1, or —OB-.—Q, is an A- in F. 
Therefore, if A has no A-! in F and B has a B-! in F, (4) cannot 
be solved by an Q in F having a reciprocal in F. In this statement 
A and B can obviously be exchanged. 

In some cases a solution of (3) can be obtained by successive 
approximations as shown in 

I. Put A = D+A®, B= 0+ B®, C= T+, where d, = ay,, 
nn» and suppose that (i) \d,—@,| >d > 0, (ii) |C| 
exists and |A|+-|B) = e < d, where the bars denote a regular bound. 
Then 


G v Pgs Yn _ c 


X=A+> x (5) 
s=1 
is a solution of AX—XB=(C, (6) 
where A, X®, X®,... are successively determined by the equations 
(D—O)A = T ) 
DxX® ~x®%@ = CO+A BO —AMA 
DX® cae X20) —_ XM BO AOXD 


DXG6+)— X69 — X® BO AMX) 


r" 


[P+ /o| 

— de * 

Proor. (7), gives 5, = y,/(d,—8®,). From (7), 

+8, Bo —a™ dj. |C|+ [Ale 
d,,—9,. : d 3 


7 


From the general equation (7), 
3) Af 0) »(s 
> xh) OY — > ay ay 
? 
d, Yes: 6. 


1 (10) | 8 s+1 
merci rag o 


Moreover, |X| < 


|X| i 


rk = 


fot-g) 
enk i. 


d d 


Therefore as 
1 ae, a | aad ee e 
iY (C@)_" 4 [A] — AM < +|A|——, (10 
z1<8 a a” 2, d—< "des (10) 
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and from the first equation 
Al = |Tld-. 


We may add that X will have a reciprocal in the field when, for 


instance, d(\C|+8'e) < y(d—e), (11) 
where y is the lower bound of |y,,|, 5’ is the upper bound of |6,,|. 

If = 0 this method gives A = 0; if C = 0 (the homogeneous case), 
the method leads to the trivial solution X = 0, for A= 0. How- 
ever, if © — D, the first equation is satisfied for every A. In this case, 
for the transformation of matrices into diagonal matrices, we may 
make use of the following remark, which I owe to Dr. R. G. Cooke. 


Sree ; : 
If X is the solution of AX—-XD=—C. 


and if C, -1X, X are associative, then 
(A—C.-1X)X—XD = 0, 

and thus X transforms A—C'.-LX into D. Some examples also con- 
structed by Dr. Cooke tend to show the efficiency of this procedure. 

Our method fails, however, for 

AX—XA=1 (12) 

which, in particular cases, can be solved step by step. For example, 
if the non-vanishing terms of A are a, ,,,,, we have 

"oer ee 5”, On n+1%n41,k—2n k-1 %-1k oe (k > 1). 


These equations leave the first row of X arbitrary; the other terms 


of the kth column are 

_  — M1,k-1 %-1,k _  _ %2,k-1 k-1k X-2,k-1 Uk-1,k 

a 2k » >. a 3k g¢es a k Lk — > 
Ayo Ao3 Ay-2,h-1 
. . . . k . _™~ ] 
Lpp = Lp_yp-1 woe = Hyy, Cpeiik ip =O (> 1). 
AKk+1 

These relations determine the kth column in terms of the elements 
of the (k—1)th column. Finally, x,,; = 1/a,., 2,,, = 0 (n > 2). 

To show the peculiar nature of (12) we first make the obvious 
remark that (12) has no diagonal solution since a,,;,8,—5,, 4, = 5% 
is contradictory for n = k. It follows that (12) has no solution which 
is the transform of a diagonal matrix, since on multiplying 

AQ“ DQ—Q-1DQA = 1 
by Q on the left and by Q-! on the right, we obtain 
QAQ-1D— DQAQ- = 1, 
so that D would be the solution of (12) where A is replaced by its 


transform. 
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Another less obvious property of (12) is the following. If B satisfies 
(12), B” satisfies ABn— BrA = nBr-1 (13) 


provided A and B belong to the same associative field F, with a 
regular bound, say. Now, take matrices A, from F such that 


lim v,/|A,| = a-, |B| <a, 


and consider the matrix function 
F(X) = > A,X’. 
v=0 


It follows from (13) that, if A, are commutable with A, 
AF(B)—F(B)A = F'(B), 


where F’(X) is the formal derivative of F(X), i.e. 
F(X) = yA, X" 
Similarly, if f(X) = 3 0. (X—By, we have 
f' (X) = S v0(X—By, 


v=1 
from which it follows that 
C, = f(B)jv'. (17) 
Therefore, if we expand F(X) = > A, X” about B, i.e. if we replace 
X by (X—B)-++ B, use the binomial theorem, and rearrange the series 
in powers of X — B, we obtain 
F(X) = 5 0,(X—By. (18) 
This process implies that we restrict ourselves to matrices X which 
are commutable with B, and for which the resulting series (18) is 
absolutely convergent at least if |X—B| < a—|B). 
If |B| < 2a, ie. if X = 0 is in the circle of convergence of (18), we 
can re-transform (18) into (14) and thus obtain 
A, = C,—Cv1C,,, B+ Cyt?C, ,.B—.... (19) 
This shows that, if C, = 0 for every v, A, = 0 for every v. Now, 
if B satisfies (12), every absolutely convergent power series with 
matrix coefficients commutable with A satisfies (15). Thus also 
F'(B), F"(B), ete., satisfy it. Therefore, if F(B) = 0, we have 
F’(B) = 0, F"(B) = 0, ete., ie. C, = 0 for every v, and thus also 
A, = 0, ie. B is not a root of an equation like F(B) = 0. 
Analytic continuation leads to some obvious extensions of this 
result. A consequence of our result is that (12), where A is a finite 
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matrix, has no finite matrix solution. In fact, every finite matrix 


satisfies the so-called rank equation. 
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1. THE present note contains the proof of the following: 


THEOREM. Let K(t) (0 << t < o) be a real, piece-wise continuous 
and not identically vanishing function with the properties 
[Kem dt << @ (every @ < 1), 
2 K(t) >0. 
We set K, .,(t) = iK, dt, K, = K. Furthermore we suppose that 


(c) K,(0) = 3. 
The integral equation o 
f(r) = | K(\r—t\)f at 

0 
possesses then one and only one solution satisfying 
f(r) 


> 


’ T > OO. 
Its value at r = 0 is 
f(0) — +4/{2K,(0)}. (3) 
For the validity of (3), (c) is the most essential part of the hypo- 
thesis. No generalization could be found in the case 2K,(0) ¥ 1. 


The particular case K(t) = } f eA dA is of interest in modern 
1 


astrophysics: (c) is easily seen to be fulfilled. The relation (3), 
f(0) = 1/v3, yields, then, the ratio of the boundary temperature to 
the effective temperature for a star in radiative equilibrium,f 
( 7) _ we 
4 ee Be 
Under the hypothesis (a), a theory of (1) has been given by N. 
Wiener and the author.{ The discussion of (1) is governed by the 


+ M. Bronstein, Zeits. fiir Physik, 59 (1929), 144; E. Hopf, M.N. R.A.S. 90 


(1930), 287. 
{ N. Wiener und E. Hopf, ‘Uber eine Klasse singulirer Integralgleichun- 
gen’: Berlin Sitzuwngsber. der Akad., 1931. 
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characteristic equation 
x(z) = [ K(t)(e*—e-*) dt = 1. (4) 
0 
Here «(z) = «(a+7y) is regular in the strip |2| < 1 and tends to 
zero as |y| > 00, the convergence being uniform in every partial strip 
2) = 2 =< 1. 
2. General facts concerning (1) 
In order to determine all solutions of (1) with the property 
- o(e%) (0<6< 1), (5) 
we proceed as follows. Let 2, 2s, ..., 2, be the roots of (4) lying in 


the strip |ja| < 6. We set 
P(z) = II e—2,) 


1 
(z? —_ 
P(z 


The function between the brackets does not vanish in |x! < @ and 


and h(z) = log}{1—x«(z)} (6) 


is positive for y = 0. The logarithm is chosen to be real for real z. 
Then h(z)—>0 for (y)+>0, |x| < @; |h(x+1y)| is a quadratically 
summable function of y. 
All solutions of (1) satisfying (5) are contained in the formulae 
O--ix 


(7) 


O+ia 


© "2 
6-—ix 


where Q,,_, is an arbitrary polynomial of the indicated degree. The 
integral (7) is to be understood in the sense of the Plancherel theory 
of Fourier integrals. Furthermore, each f(r) given by (7), (8) vanishes 
for + < 0 and is continuous for 7 > 0; (7) can be inverted, and so 


F(z) = [ f(ter dt (w < 0). (9) 


( 


From (7), (8) the asymptotic form of the solutions can be ob- 
tained,t for 7 large, 
fir) = ¥ Q*(x)e*"+-0(e), 
+ By transferring the integration abscissa in (7) to the right as far as 0 and 
by taking account of the residues of F'(z)¢ 
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where z* is a root of P(z), and where Q* is a certain polynomial of 
degree less than the multiplicity of z*. 


3. Proof of the theorem 


We have «(0) = 2K,(0): (c) means therefore that z = 0 is a root 
of (4). Since «(z) is an even function, we have «’(0) = 0. Further- 
more, Fs 

«"(0) = 2 [ #K(t) dt = 4K,(0) > 0. (10) 
0 
Therefore z = 0 is a double root of (4), and, according to (b), no 
other root can lie on the imaginary axis x = 0. The roots of (4) can 
thus be arranged in the form 


~ 


Zy 900+) Sy —3, 0, O, —Zq5.--) —Z yao 
where —6< 2, = Ri(z,) < 0. (11) 
n—1 
On writing P*(z) = JT] (z—z,), (12) 
1 

we have P(z) = (—1)"—122P*(z)P*(—z). (13) 
If, now, (1) has a solution which is approximately linear for 7 large, 

we infer from (9) that F(z) is regular for x < 0. This regularity can 

be obtained only by setting 


Q,,-1(2) = cP*(z) 
in (8). We have thus 5 
” ne @—" exp) 1 ff WO 
F(z) = (—1) toot e| ont [ “ae ac (14) 


6-ix 


Conversely, we prove that (7), (14) actually yield an asymptotically 
linear solution of (1). F(z)e-7* can have poles only at z = 0, —2,,..., 

z,,. Only 7 and bounded terms can therefore occur in the asymp- 
totic expression for f(r). For ¢ #4 0, however, 7 cannot be missing, 
since a bounded f(r) would, according to (9), imply 


|zF(z)| < const. |z/n!| (x < 0), 


a contradiction to (14). 
For y = 0, x < 0, we infer from (9) that 


la| F(x) = [ fi(t/\x\)e~ dt, 


0 


thus yielding || F(a) > f(0) (a > —oo0). 
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On the other hand, from (14), 
|u| F(a) > (x > —oo), 
thus yielding c= f(0). (15) 


Now we investigate F(z) at z= 0. If (2) is satisfied, we know 
that f(r) = 7+ O(1). According to (9) this implies, for z < 0, y = 0, 
f l 
P(x) = ae O(1)e™ dt = —+ (=), 
ge CS x? x 
0 . 


thus yielding lim 22 F(x) = 1. 


x—>—0 


On the other hand, from (14), 


c—>0 Pr 27 4 E 
O—ixw 


O+ix 

or a! |, ee ] " h(€) * 

lim z?F(z) = = e%?|—3 : | ie (17) 
The integral in (17) may be evaluated as follows. The path of 
integration, R(é) = 0, may be transferred to the left and deformed 
into another broken line, which consists of the two parts (— 100, —ip), 
(ip,i00) of the imaginary axis and of the right half of the circle 
£| =p. It is now essential that the integrals along the first parts 
cancel each other, for h(€)/€ is an odd function. Hence 


0+ ix i 


[ Mae = + F aie a 


. 


This tends to 7zih(0) as p > 0. Therefore (17) becomes 
“ re Cc , 
sap nia Peo)o 

thus yielding, according to (15) and (16), 

f(0) see P*(0)ei*O, 

From (6), (10), and (13), however, 

h(0) = log $x”(0)—2 log P*(0), 


the logarithm being real, and thus (18) gives the result (3) 
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By L. CARLITZ (Cambridge) 
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Introduction 
IN a paper in this journal,} bearing the same title as the present one, 


I considered the sum 
a(;)...0(n,), (1) 


N=N,+...+Ny, 
taken over all the partitions of n into v positive integers, a(n) being 
a ‘factorable’ arithmetic function satisfying certain order conditions. 
By applying the Hardy-Littlewood method, an asymptotic expres- 
sion was found for (1), at least for vy > 3. The present paper deals 
with the somewhat more general sum 


> %4(1).--4,(m,), (2) 
for all values of v > 2, where the functions «;(n) are similar to a(n) 
above, that is, they satisfy the same general conditions that are 
satisfied by a(n). The method used is now quite elementary. It was 
suggested by Estermann’s{ proof of a theorem due to Evelyn and 
Linfoot§—a theorem included as a special case of the theorem to be 
proved about (2) (for vy = 2). We find that the sum in (2) is 
Bn’ )(n) + O(n’-*+), (3) 
for all positive «; where B is independent of n, but depends on the 
functions «;,(v); a(n) depends on the a«;,(n), and is a function of 
exactly the same kind. While the functional dependence of a on 
v; is fairly simple, the expression of a”) explicitly in terms of «; is 
quite complicated. Accordingly, we postpone the precise statement 
of the theorem to be proved until we have made the necessary 
definitions of certain sets of functions derived from the «; (see § 1). 
Instead of considering directly the sum (2), we find it convenient 
to consider the more general sum 


mito (,)...miva,(n,), (4) 


N=N,+...+N, 
where the h; are non-negative integers. The generalized sum (4) is 
of some interest in itself; the point of the generalization is, however, 


+ Quart. J. of Math. (Oxford), 2 (1931), 97-106. Cited as I. 
t Journal London Math. Soc. 6 (1931), 37-40. 
§ Journal fiir Math. 164 (1931), 131-40. 

3695.3 , 
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that once the case v = 2 has been treated, we are able to pass 
immediately from the theorem for v to the corresponding theorem 
for v+-1: furthermore, the presence of the h; adds very little com- 
plication to the treatment of the case v = 2. 

The case v = 2 of either (2) or (4) naturally suggests the conjugate 

sum ¥ a, (n)a(n+h). (5) 
noe 
This sum can be treated by the method employed for the sum (4). 
If k be fixed, then, for x > o0, the sum in (5) is 
Baol?(k) + O(xt+*), (6) 
where B is free of both k and x, but depends on «a, and a,, and a 
is identical with the corresponding function in (3) (v = 2). 

We next give some simple applications of the theorems concerning 
the sums (4) and (5). Finally, we consider briefly the hitherto excluded 
case v = | of the sum (3): in place of an asymptotic formula we now 
get an identity. This identity suggests several others that are rather 
curious, and seem to be of some interest. It should be remarked that 
for special values of a(n), the identities obtained reduce to identities 
discussed by Evelyn and Linfoot (see § 6 for reference). 

1. Definitions 
Let a(n) be an arithmetic function which is 
(i) factorable: that is, 
a(mn) = a(m)a(n) for (m,n) = 1; 

(ii) ‘linear’: that is, 

aft} ==-3. a(n) = O(n*) (e > 0), 

a(p) = 14+ O(p-**),t 

where p, as throughout this paper, denotes a prime. 

Let B(n) be the so-called derivative of a(n): 


B(n) = 2, H(d)a(5), 


where, as usual, u(d) is the Mébius function. It is clear that 8(n) has 
the following properties: 
(i) B(n) is factorable; 
(iii) 6) = 1, B(n) = O(né), 
Bp) = Op") (e > 0). 


t The } in the exponent might be replaced by any positive quantity ; how- 
ever, for simplicity, we limit ourselves to the case that gives the best error 
term in our theorems. 
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We shall say for brevity that any arithmetic function satisfying (i) 
and (ii) is an alpha-function; any function satisfying (i) and (iii) is 
a beta-function. 


Lemma. If B(m) is any beta-function, then 
P B(m) = O(at**). (7) 


Define the function @(m) thus: 
(1) = 1, 
6(m) = 1 if for p|m, p?|m, 
6(m) = 0 otherwise. 
Then we may verify that, for s > }, 


i a) 


O(m) __ £(28)¢(3s) 
Ps = a = (i) a — 


whence it is easy to show that 


Y A(m) = O(z'); (8) 


and indeed much more may be proved by elementary methods. 
Now it is evident that any integer m may be written in the form 


m = gh, where p*(q) = 1, 0(h) = 1, (¢,h) = 1. 
Accordingly, 


5, om)| =| _oHm)ea) 


m=1 


Max 
(q,h)=1 


= =i \B(q)|9(h) |B(A) | 
~ og, 3, 000 


h<a/q 


- Olz™ > ee/a)'), by (8), 


= Olt ¥ 4} = O(at+€), 


This proves the lemma. 
Using (7), it is clear that the series 


b(m) = S om) 


is absolutely convergent. We put 
b(m) 
T(m) = ol)’ 


T2 
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We shall now show that 7(m) is a beta-function. Evidently 
b(m) = Ty {apr 4 PP) + |, 
p \ Pp J 
where m = [| p*, the product being taken over all p (so that 
e = e(p) = 0 for pim): it will be convenient to use this notation 
throughout. Then, by (10), 


ppt) + OP) 


(m) =| | OF (11) 


pim 1 "EP cen 
Pp 
(We assume, of course, that A(m) is such that no denominator in 
(11) vanishes; that is to say, b(1) 4 0.) From (11) it is obvious that 
t(m) is factorable. It is also fairly clear that 7(m) has the other 
necessary property, (iii) above. For the sequel, the relation between 
the beta-functions, B(m), 7(m), is important, and we shall write, as 


an abbreviation for (11), —) (12) 





As for the inverse B = Q-17, it is easy to see that, from (8) or 


(9), we get = scfe\r(ona 
B(m) = = where t(m) = > u(s)2 (me) (13) 
s=1 


s 

e+1 

T(p*) m. ) 
so that (m) = ee aS (14) 

Pm): I] 1 —7(P) 

Ps 
We remark that b(1jt1) = 1. (15) 
It is clear from the definition of a beta-function that the product 
of several beta’s is itself a beta; also m~18(m) is a beta whenever 

B(m) is one; therefore, plainly 


7) = m'-"7,(m)...7,(m) (16) 


pim 


is a beta, where, of course, 
7 = OB; for i = 1....,v. 
Then, as we have seen above, the function f™, defined by 
p™ = Q-1,™ 
is a beta, and finally, therefore, 
v(m) = + B(m) 


dim 


is an alpha-function. 
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If we write b[f] for b(1), defined by (9), and ¢[7] for ¢(1), defined 
by (13), or more briefly, 

b= O68), @M= qr), (18) 

then we may state our first principal theorem in the following form. 

THEOREM 1.f If «,...,a, be v alpha-functions; h,,...,h, be non- 
negative integers; then, for all vy > 2, « > 0, 

Rn) = —  _ mfay(m,)...nhra,(n,) 


N=MN+...+n, 
= H,b,...b, t natty y-1g)(n) + O(na tthe tv-4t+4), (19) 
where a(n), b;, are defined by (16), (17), and (18); and 
a Gate-thto—D! 
? h,!...h,! 
We proceed now to prove the theorem first for v = 2, then, by an 
easy induction, for all vy > 2 





2. Proof of Theorem 1 for v — 2. 


2.1. From the definition of B;(m) we see that 
R(n) = R,(n) = Be _,Pal €1)Bo(€2)(d &)""(de eg)". 


d,e,+ 


Let R'(n) = > ‘Bi (c,)Ba(es) (d; €,)"(dye,)". 


n=d,e,+d, 
1, €, <n 


Then — | R(n)—R(n)| <n > +S -)(Bx(ex)B2(e2)|- 


1>nels é>n* 


But > IBA (€,)Bo(€s)| = wt IBi(e:)| > Palea)| 


n=d,e,+d,¢. d,e,=n—d, 
eé,>ns Ses 


= O(n) ¥ |Bx(e)| = 
e>n*8 
[for > |8(d)|, where d|m, is an alpha-function of m, and therefore 


= O(m*)] ~ O(n") S 1B, (e)| 


e>n?ls 


= O(1 ni+€), 


by using (7). Therefore, 
R(n)— R'(n) = O(no+te+t+e), 
2.2. Now, from (20), we may write 
R'(n) => > By (€1)Bo(€2)(dy €1)"" (dy en)” 


n=d,e, +d e, 
(€1,€:)=M; 41, @<nil3 


= ¥ Smh, ny 2 
in 


+ It is easily verified that, for a, =... =a,,v > 3, h, = 0, this theorem 
reduces to the result of I; the present error term is slightly better. 
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where in = = (4, fx)" (dof2)"By(mf,)Bo(mf2) 
tims rf Sie Se 
= DL SPP mf) B(mf2) SY dade. (23) 
Y = 1 ; n/m=d,f,+dsfe 
1 n*)3 m- 


To evaluate the inner sum in (23), note that 
d, = d; (modf,) (0<d, < fy), 
d, = d (mod f,) (0<d,<f,), 
so that the sum in question becomes 


fpf > ‘ (8, +9;)""(5,+65)", 


n’=8,+8, 


where a” we oe +0(1) (0 < 0,, 6, < 1). 


mite 


But, by the Euler-Maclaurin sum-formula or otherwise, (24) reduces 


” h hs | h, th.! 2° nN Ii+he+1 n hy the) 
i : hs \( } 7 — O -- 
thet DW! \mps fe mf, fs J 
oe hi+he ie ee a! n hi+he : ‘ ) 
A | “hi a (=) fits if 
Substituting from the above, (23) becomes 


S' =H, (2) ho+1 S By(mf,)Bo(mfe) 
m fle 


m 
(fife)=1 
Sofe<n?*m" 


hi+he 
+0(*] M3 |B, (mf)Bo(mf.) | 
Sutecn**m-* 


But > |B(fm)| = O(mntm-1)i+*« = O(ni+*), 


S<v?Bm7 


hy +he+1 
sin j ss 
Thus S, = H, ( ) Sr 4 ( 
m j 


gn B,(mf,)Bo( (mfe) 
“oS fi ra 
(fi,f2)=1 


Sir ficnv?m- 


B,(mf;)B. (mfz) 
fife 


where 
2.3. We now let 


fis fe=1 
(fi,f2)=1 


then <i By we > By el 


fi>n*m- te 


> |B ( pul > pal 


f,=1 I f,>n*3m> 











ON A PROBLEM IN ADDITIVE ARITHMETIC 


But S Be ~ br*(m) = O(m*) = O(n®), 


where 7* = Q|8|; on the other hand 
] m +e t+e 
T° ee ae 
f>n?m> f>n**m- 
— Of{m*+«(ntm-!)*+*} + O(n-t+€m1+*) 
= O(n-*+€m1+«). 
By (26), therefore, 
S,,—Sm = O(n-**«m}+**). (27) 
2.4. At this point we pause to combine the results of 2.1, 2.2, and 
2.3. By (21), 
R(n) = Ri(n) + O(ns+e4+4) 
= ¥ mates) +O(nute+tte), by (22), 
pe 
a hi+hz+1 <m hi the+t+e 5 
=H,n y 3 m+ O(n ), by (25), 


_— H, ginthet1 > Sim O(nlothett > n-§+€m*)+ 
m 


min 
m\n 


+O(nlrtheti+e), by (27), 
S ae 
a Hy, qinthe +1 p 3 a O(nlathett+ ). (28) 
mn 
Clearly, then, by the remarks made at the end of § 1—in particular, 
by (17)—in order to complete the proof of our theorem (v = 2), we 
have merely to show that 
S,, = 6, b,t?mB(m). (29) 
This we now verify. 
2.5. By the definition of S,, in 2.3, we have 


‘. Sins sas e B,(mSf,)Bo(maf.) 


P 2 
d=1 8 §=1 (fysf2)=1 Shi fe 


S B,(mf,)B2(mf >) 
Suvf,=1 fifs 

= b,7(m)b,72(m), (30) 

as follows from formula (10). But inverting (30), we have immediately 


ao 


(8) 
S,, = 6,6, —-74(md)r_(md) 
2's 
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= b,b.m re (9) (md) 


= 6b, b,t? mB (m), 
where 
72m) : LE Sead in B®) = -17(2), {2) — [7] 
m 
[see (16) and (18)|. We have then proved (29), thus completing the 


proof of the theorem in the case v = 2. 


3. Proof of Theorem 1 for v arbitrary 
3.1. Let us assume now that Theorem | is true up to and including 
the value v. Then we have 


Ryas(m) = XB le)aysal ff" 


n=e 
= H.b,..b,1 D> elite thy ty-lyM(e) fletig, (f+ 
ae > to thotitv—t+e), (31) 
the error term being estimated exactly as in the case of the sum 
in (24). 
Since «(e) is an alpha-function, the sum in the right-hand member 
of (31) is an instance of (19) with v = 2. We write 
h,+...th,+v—1, | a”), Oy445 
in place of hy, A,, Oy, Xs 
respectively. Then the sum in em is 
m4 (hy +... th, +v—1)!h,,,! 6B BEB, aa Mi 7*] 


4. iat 


n=e+f 
SS *(n)+ O(nla tts ritv—$+e) 
where ofp” | = Bote aes by (15); 
eal ~ 4)? . ? 
7*(m) = T(m)ty41(m) _ 74(m)...7,(M)r, 44(m) 
m m” E 
so that 7* = 7th), {r*] = te 
finally *(m) = > B*(d), p* = Q-17*, 
dim 
so that B* = Be and a* = a’t, 


Substituting into (31), it is clear that the induction is complete. 


3.2. We wish here to make two remarks about Theorem 1. In the 
first place, we note that the more general case in which the h; are 
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merely real and each > —1 is by no means ruled out. The only 
places in which the proof is affected are those in which the Euler- 
Maclaurin formula is used. If, then, we make one obvious modifica- 
tion in the statement of the theorem, namely, 


yw — Pat+D)...Ph,+1) 
"Py +...+Fh,+y) ’ 
we may still assert the truth of the theorem, at any rate with the 
less exact error term 
O(n +-+%v+¥-1) for some c > 0. 

Secondly, by modifying the proof in § 2 somewhat, it is not difficult 
to see that a direct proof of the case v arbitrary may be given. 
Xestricting ourselves, for simplicity, to the case in which all the h; 
are zero, we see that the only difficulty is the evaluation of the inner 
sum in the general form of (23), or, what amounts to the same thing, 
the determination of the number of solutions in positive integers x; of 
nN = 4,2,+...+4,2,, 

where (4,,...,a,) = 1; the number of solutions may be shown to be 
. nY-1 n-2 
a = ee o(-— 2 a4). 
[ hope to show elsewhere in detail how this formula may be employed 
in a problem similar to that with which this paper is concerned. 


4. The conjugate sum for v = 2 


We now consider the sum (5)—in slightly more general form 
involving the h’s. 
THEOREM 2. If k is a fixed positive integer; h,, hz are non-negative 
integers; a4, % are arbitrary alpha-functions, then 
R(x) = ¥ noay(n)(n-+k)hrxg(n+k) 
nor 
eis b, bt (k)+ O( hath ++) 
= —___ a ths , 
h,t+h,t+1'°* 
where b,, by, t, a are defined as in the case v = 2 of Theorem 1. 
4.1. As in § 2.1, we begin with 


Rie) = Y — Byleq)Bol€a)(dy e1)"(d ee)". 


dye,<r 
de, +k=dye, 


Let R(x) = > By (€1)Bo(€2)(dy €1)""(dp ep)". 


1€, + K=A:C, 
dye, <2; €;,€, <2" 
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Then | R(x)— R'(a)| < x(a-k)( a7 >) Br(er)B2(e2)I. 


But >, |Pa(ex)Bo(e2)| on ‘yi Pr(es)| | : > Peles)! 


= Ofe*) ¥ |B(e)|= = O(at+*); 


and similarly 
S, !Pa(ex)B2(¢)| = 1 2 - B,(¢2) > B,(e,)| 


| 
d,e,—k=d,e, 
ey > > 2 '3 


Therefore, by (32) and (33), 
R(x) —R' (x) = O(alattatt+e), 
4.2. We now write 
R'(x) = ¥ Sj, mb, 


where Sn = - (4, f1)"" (de f2)"B (mf, )Bo(mf)- 
k/m=dzf,—d,f, 
(fupfe)=1 
AA, <23 fir f2<e*m-* 
Then proceeding exactly as at the corresponding point in § 2.2, we 
find that 1 ( 


«3a 


A 


sate ml 
m 


where Ss” = > By(mfs)Bolmfa) 
(fisf2)=1 hfe 
Sishe<x??m-* 

4.3. If we put | a By(mfs )Bo(mfa) 


: fife 


(fi, f2)=1 


v hit+he+1 ’ e 
Si, O(ahrthertre), (omy 


then, precisely as in 2.3, 
S,,—S), = O(a-***), 
4.4. We now combine the results of 4.1, 4.2, 4.3. 
R(x) = R’(x)+O(ah+at+i+e) by (34), 
= Y mates) + O(ginthattte) by (35), 


milk 
ginths +1 


‘ ‘mn : 
= Sm 1 O(girtheti+e) by (36) 
| T A ? 
hy+h.+1 = = 
aliy+he+l S 3 
iz - : : m4 O(aitati+e) by (37). 
“rT tet m\k 1 
4.5. Evidently, in order to complete the proof of the theorem, it 
is only necessary to show that 
S,, = b,b,t?mB(m); 











ON A PROBLEM IN ADDITIVE ARITHMETIC 283 


but this is simply (29), and the quantities involved in the two equa- 
tions have exactly the same significance. Hence the result of 2.5 
may be applied, and the proof is completed. 


5. Applications 


5.1. Let the arithmetic function p,(n), where s(> 1) is a fixed 


integer, be defined by 
p(n) =1 for p*jn, 


w(n) = 0 for p*|n; 
it is easy to verify that u,(m) is an alpha-function. We apply the 
definitions of § 1 to calculate the corresponding f and 7. 
(i) B(p*) = —1 for ¢= 6, 
B(p*) = 0 for l1<e<sore>s. 
(ii) t(p*) = 0 for e>s, 


— 
7(p*) = - P for l<ec<es. 
pe—l 
Let now §8},..., 8, be a set of vy integers each greater than 1; assume 
them so numbered that s, is not greater than any of the others. De- 
fining the set of alpha-functions 


aj(m) = Hs,(n), 


we shall now calculate 7”, 8, and a. 


(iii) mp’) = 0 for e>8,, 
(—1)"p° 
7)(pe) = ___"2______ for 
a (p**—1)...(p*—1) 
(iv) Bp?) = 0 for e~s, (e>9), 
(-Yr 
(p—1)...(p**—1]) 
v)(me) — 
B(p*) pea (—ly 
(p*\—1)...(p**—1) 
wp’) = 0 for loe<s,, 
Bes TAN 
ape) = an for e>8,. 
~ (ph—1)...(p—1) 
We may now state the theorem. 
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THEOREM 3.7 If hy,..., h, are non-negative integers; 8,,..., 8, integers 
each > s, > 1; then the sum 


> n..n'v, where p*ifn;, v > 2, 
N=Ny +. Ny 


. H, Cy(n)nlo—+ »+v 1 O(niat—the tv—§ +€), 
,!...h,! 
(h,+ eee th,+v—1)!’ 


| 
| T'- ~ (p—1)...(p*—1)J’ 


oe a ee 
(p*—1)...(p*2—1)} | (p*—1).. (p*—1) 


pin ~ 


where B= 


and 


n(n”) 


5.2. We now consider a more general problem. Let s(p) be a func- 
tion of the prime p, whose values are integers not less than 1 always, 
and indeed not less than 2 except for at most a finite number of 
primes. Define the factorable function g(n) by 

g(1) = 
g(p’) = 1 for e< s(p), 
g(p?) = 0 for e>s(p); 


it is easy to see that g(n) is an alpha-function. The corresponding 
B and 7 are defined as follows: 

B(p*) = —1 for e= s(p), 

B(p*) = 0 for 1<e< s8(p) ore >s(p); 


7(p*) = 0 for e > s(p), 
€ < 8(p). 


Let now 8,(p),..., 8,(p) be a set of v functions of the prime p, 
satisfying the conditions satisfied by s(p) above. With each s,(p) we 
associate an alpha-function g;: we have already calculated the corre- 
sponding f; and 7,; we shall now calculate 7™, B™, a, 

+ This theorem in the case h, = 0, s; = ... = 8,, is due to Evelyn and Lin- 


foot, Math. Zeitschrift, 30 (1929), 433-48; see also I, p. 105. The case v = 2, 


h = 0, is due to A. Page, Journal London Math. Soc. 6 (1932): from this 
paper it will be clear how the error term in Theorem 3 may be improved. 
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(iii) 7+(p*)=0 for e> s(p), 
(—1)*p* 
(p*™—1)...(ps*)—1) 


where s(p) = mins,(p); 
1<i<pv 





7(pe) = for 1<e< a(p), 


(iv) Bp?) = 0 for e~s(p), e>9, 


(—1)"p° 
B® (p*) ae (gn 1 - —_ 1) for 


(ps) —1)...(p*)—1) 








(v) a”(p) = 0 for 1<e< s(p), 
(—1 fig) 
(p* —1)...(p“)—1) 
(—tF 
(ps) — 1)...(p“)— ] ) 


+ 








of) pe) = for e > s(p). 





THEOREM 4. If h,,..., h, are non-negative integers; 8,(p),..-, 8,(p) 
functions of the prime p such that (for each i) s;(p) ts integral and > 1 
for all p, and s,(p) > 2 for all p with at most a finite number of excep- 
tions; then the sum 


> nh.ny, where pn; v > 2, 


N=Ny +... Ny 
a H, Cn(n)nlnt~+hetv-1 4 O(qlat~tty—t+e), 
where H,, is as in the previous theorem, 


Pi, (—1)" — 
; Tp ayy | 10 pP ))..(1 p MP), 


p p 





and 


(—17“gi—8 (—1y a 
1— 1— - 
7(n) P| es 


py Pin 


8(p) being defined by (38). 








5.3. Let o,(m) denote the sum of the sth power of the divisors of 
n; evidently, for s < —1/2, o, is an alpha-function. For s merely 
negative, the corresponding function is one for which Theorem 1 still 
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applies, but with inferior error term. It is easy to verify that 
Theorem 1 implies the following result (see the first remark in 
§ 3.2). 

THEOREM 5. Let 8,,..., 8,; h,,..., h, denote (real) positive constants; 
then, for all v > 2, 


> = nhbo_,, (m,)...np”-10_, (n,) 


N=N+...+Ny 
~ (8+ 1)...0(8,+ 1) P(h,)... (hy) niit~—+thy-lg_ (mn) 
C(s+1) D(h,+...+h,) aah 
where 8s = 8,+...+8,+r—l. 

5.4. Theorem 5 is capable of a simple generalization: the generaliza- 
tion has the same relationship to the theorem that Theorem 4 has 
to Theorem 3. 

Let s(p) be a function of the prime p such that 

s(p) >¢ > 0,4 
where c is independent of p. Then we define an alpha-function (in 
the more general sense: see § 5.3) by 
oa’ (p*) —_ 1+ pM)... n—es(P) 

THEOREM 6. Let hy,..., h, denote positive constants; let s,(p),..., 8,(p) 

denote functions of the prime p such that 
s(p) >¢ > 0,F 
where c is an absolute constant; let o';(n) be a set of v factorable functions 


of n defined by o;(p*) cnet l+p s(P) vee pep); 


then the sum 


nir—..myr2a;(m)...0,(Ny) ~ 


D'(h,)...P'(h,) 7 na 
Re as | te v B hyt+...+h,—1 : 
r@,+..+h,) a(n), 


. > - ts l = p 8(p)-1 
where B= {| ET — ae 
(l—p )...(l—p 


p 
o’(n) is factorable, 


o'(p*) = 1+ p3)+...4+-p-e), 


and 8(p) = 8,(p)+...+8,(p)+v—1. 


+ It is sufficient to assume 
loglog p 
s(p) >e ES , where c > 0. 
log p 
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5.5. If d(n) is the Euler ¢-function, then obviously ¢(n)/n is an 
alpha-function. More generally, if ¢,(m) is Jordan’s generalization of 


the ¢-function: $,(n) = n* TI (1—p-*), 
pin 


then n-*¢,(n) is an alpha. Hence the following theorem: 


THEOREM 7.7 If hy,..., h,, 8,..., 8, denote positive constants, then 
the sum 


nh—8.-1, a) hy t—1d (m,).. ds, (n,) 


D\(hy).--P(hy) ate thy— 
“he. ee 


where 





BE(s,+1).. £(s,+1) =|) Tt — Gah 


and y(n) is an alpha-function defined by 
ea ee ee 
(P= 1)... (A= 1)..(pFA—T) 


pin 








It is clear how this theorem may be generalized by replacing each 
s by a function s(p) that is greater than some positive constant: in 
view of the close analogy with previous theorems it is scarcely neces- 
sary to state the generalized theorem. 


5.6. So far we have given applications of Theorem 1 only. How- 
ever, comparison of Theorems 1 and 2 shows that for any application 
whatever of Theorem 1, the corresponding application of Theorem 2 
may be written down immediately and without any further calcula- 
tions. Therefore, we shall not take the space to enunciate the 
theorems conjugate to the case v = 2 of Theorems 3-7. 


6. Identities arising from the case v = 1. 

We have had no occasion in this paper to introduce the ‘singular 
series’ for the sum (2). It appears likely from the result of I that 
this should be 

g a = 7k) ae . 
un) = > k >= > ae eS 
1 p(k) k=1 p(k) 
where p = p(k) runs through the primitive kth roots of unity. Indeed, 
following I, § 4, it is easy to prove that 
S,(n) = Ma(n), 
+ For v = 2, see A. E. Ingham, Journal London Math. Soc. 2 (1927), 208. 
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so that the conclusion of Theorem 1 may be written (for h; = 0) in 


the form 


R,(n) ~ at S,(n)n’-} 


We remark that in the case v = 1 this asymptotic relation reduces 
to an identity. It is only necessary to show that 


M—t, and a” = a. 
The first equation is obvious from the definition of ; as for the 


second, since 7 = 7,, it is evident that BY = f,, and therefore 
x) = «,. Dropping subscripts, we put this result in the form 


a(n) = bS(n) = b - 7(h) ps 9. (39) 
k k - 
K=1 p(k) 
There is no difficulty in showing that the series converges. 
Let now f(x) be the generating function 
fle) =X a(nye"— (|| <0); 


then, formally at least, from (39) we get 


‘) 


ei eel 
f(x) = bx , 


k p(k) ts 


This can be proved directly in the following way: 


¥ a(n) x” = Fx" ¥ Bd) 


n=1 n=1 d/n 


. B(d)at 


l — yt 
d=1 


iz => a> t(ds) 


d=1 


SS r(k) S det p(s) 
= by > S > ; a ae (41) 
c=1 k=ds ‘ 


Let F(x) denote the polynomial of lowest degree whose zeros are 
the primitive kth roots of unity; to make F, unique, assume the 
term free of x to be +1; then, as is well known, 


Fy(e) = T] (1—wt wn, 


k=ds 
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It is obvious, then, that, using (41), 


which is equivalent to (40). 
We next consider the sum 
— (ke 
b, > BOS pup), where fix) = 
k=1 ~ p(k) 
By (39), this sum 


=: > a4 (n)ag(n)a" = > a(n)x”, 


n=1 n=1 


where a (= a, 49) is itself an ek We have, therefore, 


n 


p(k) 


= b, ee . (43) 
k=1 ; 
This result reduces to (40) on putting a,(”) = 
As an illustration of these formulae, let a(n) = g(n), where g(n) is 
the function defined in § 5.2. If we write g{n;s(p)} in place of g(n), 
then formulae (39), (40), (42), (43) reduce to the following direct 
generalizations of the Evelyn-Linfoott identities: 


gin; s(p)} = b > om 8( +9 | | am pol Pg (39’) 
jo) = Sotmoenper = we > oto +Y] oa DS Ow) 


= —be Sotto] Tea Hay 2) 


+ Annals of Mathematics, 32 (1931), 268-70. There are several other 
identities given by Evelyn and Linfoot that may be immediately generalized 
to the case of a general alpha-function. 

3695.3 U 
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and 


Sol: s(p)} a = =, > ako $1(P) ‘at | pi? — 5 > fale) 
1 


= bs > otkssr+3 | | a Ty Dale), (43") 


k=1 pik p(k) 
where : II (1 —p*?)), b; II (1 —p Si(P)) 
p p 


and in (43’) s(p) = mins,(p). 


1=1,2 
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By A. G. WALKER (Edinburgh) 
[Received 11 March 1932] 


1. Introduction. This paper is mainly a discussion of the first and 
second curvatures of a Riemannian sub-space. The first curvature 
is well known as the mean curvature, and has been examined by 
Eisenhart and other writers; the second curvature is introduced 
directly as a generalization of the squared torsion of a twisted curve. 
The notation employed is, in general, that used by Eisenhart,* and 
we shall quote many of his results without giving explicit references. 


2. Notation. Consider a flat space V,,, coordinates 2 (t = 1, 2,...,p). 
Let V,,, coordinates y* (a = 1, 2,...,m), be a sub-space of V,, and 
let V,,, coordinates x‘ (i = 1, 2,..., n), be a sub-space of V,,. Then 
V,, V, are given by equations of the form 

g=Hy); g*=yg{z), Oe = Te), (¢ = 1, 2,...,p). (1.1) 
As J, is a flat space, we can use the generalized vector analysis when 
referring to this space.t Let r be the position vector of a point of 
V,. Then r is a function of the y’s for points of V,,, and of the 2’s 
for points of V,. We shall use the Greek letters a, 8, y, 5, and the 
roman letters h, i, j, k for suffixes when referring to V,, and V, 
respectively. Thus we shall write f, = éf/éy*, and f; = éf/éa*. Also, 
fg and f;; will denote the second covariant derivatives of f with 
respect to V,, and V,, respectively. 

The fundamental tensors of V,,, V,, are given by 

ag = r, Tp, 9:3 = T;°T;, (1.2) 
and from the relation r; = r,y%;, we have at once g;; = d,gy%; y? 5. 

There are m—n independent normals to V,, in V,,, and p—m 
normals to V,, in V,. We shall choose these to be mutually ortho- 
gonal, writing them N,; (o = 1, 2,...,.m—mn), and N,, (r = m—n-+1, 
...,p—n). Thus we have 
No N,, =9 Nj =¢e, = +1; N,,-N, = 90, N? =e = +1, 
and N,)'T; = 9, N,,T,, = 9, (1.3) 

(o,n = 1,2,...,.m—n, 17,8=m—n+l,....p—n, += 1,2,...,m, 

e= i, 2....;m, of~yp;r #8). 

* Riemannian Geometry (1926). 


+ For an introduction to this use of the vector analysis, see the writer’s 
paper ‘On deformation of sub-spaces’: Proc. Edinburgh Math. Soc. (1932). 


U2 
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As the normals N,, are tangent to V,,, we can write 
pe . 9 
N., = £99. (o = 1, 2,...,m—n) (1.31, 
where &,,~ are the contravariant components of these normals in VJ,,,. 
The second fundamental tensors of V,, in V,, are given by 
Ooi F N, Ti aa —Fy ty — — ais 8 | 


omit N, ; Nui hag Nous N,. 


For V,, in V,, we have the above tensors, together with the tensors 


(1.4) 


pb 


Q...;;, Hysi¢ Similarly defined, and also the tensors 
Horii = N, Nive (1.41) 


From (1.4) and (1.41) we see that the fundamental equations for 


may now be written 


r ij Dd ee Qoiizi No + ¥ e,0,;;N, | 


n ) 


V, in V, 


5 rij r 


Go 


(1.5) 


a,t 
T 


= 
ik wr J 0 4 
N } Q, 7’ Pr, ae p 3 Cu Fou iN, = > Cy Porli N, 
pe ] 


and the equations for V, in V,, are found from these by omitting 


vt 


vector components normal to V,,,. 
: 7 2 aw - ol 
For V,, in V,,, we have introduced the vector operator V = gr; —. 
ox! 
We can now extend this for V, in V,,, and define the relative operator 
V — gir Om 
=e Seno 
COX 
where @,,/0x! is é/éa/ followed by an operator for annihilating vector 
components normal to V,,.* In this way, we shall extend the operator 
V2? = V-V, ete. 
Lemma. 1'o show that 


giiy* y ; =—2 > e, é, xe, B. 
o 


Let A,’ (7, = 1, 2,...,2) be the components of an orthogonal en- 
nuple in V,,. Then the components in V,, are A,,* = A,,‘y*,, and these 
n vectors, together with the (m—n) vectors €,,*, form an orthogonal 
m-uple in V,,. Hence, we have 


a*B — ¥ e,A,°A,F + > €, €,:%E,F. (2.1) 
h a 
But > en rj, *A,, B = a yP ; bY en rj, , f= gly; uP (2. 1 }) 
h h 


t 
* For a discussion of the vector operator, see the writer’s paper above 
mentioned. We shall use the same notation V for the relative operator: for 
this is the only one with which we shall be concerned. 
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Therefore gly y® ; = a8 — > e,é,,%E,P. (2.2) 

o 

If m = n+1, this becomes 

ExgP — e(arP—giiys ,yB 5) = cP, (2.3) 
say. Hence the components of the normal to V,, in V,,,, are given by 
&% + e,(c%*)H2, (2.4) 

where e, = 1, and e, is the sign of c!*. 

For a surface in V3, (2.4) becomes, with the usual notation, 

X = {1—(Ga2—2F x, x,+ Ex2)}"2, ete., (2.5) 
where x, = 0x/0u, x, = Ox/év. 

3. For a curve in V3, the vector operator reduces to V = td/ds, 
where t is the unit-tangent vector. We find that the curvature vector 
of the curve can be written 

cn = Vr, 
where n is the principal normal, and « the curvature. Generalizing 
this expression to V,, in V,,, we write 


MN = V°r = git ¥ e254; Noi: (3.1) 


We see that N so defined is the mean-curvature normal, and M is 
the mean curvature of V,, in V,,, these, then, being the generaliza- 
tions of the principal normal and the curvature of a curve. We shall 
call N the principal normal, and M the first curvature of V,, in V.,,. 
The second fundamental tensor for the normal N we shall write as 


Q;;. Hence, 


MQ,; oa b 3 €,.Q, Qos Q, wae J Qo ii5s (3.1 1) 


and 


M = &Q, M* = e} ¢,Q,); Q = g'Q;;, N?=e= +1. 


Also, for a curve, we find that the torsion is given by 
7? = —n:-V*n—(V-n)?, 
where n is the principal normal. Hence, generalizing, we shall de fine 
the second curvature T of V,, in V,, by 
T = —N-V?N—(V-N)?, (3.2) 
where N is the principal normal of V,, in V,,._ We have* 
N- VN = g@N-[Nijln = —9[NG: Ni lm (3.21) 
* The suffix m denotes that, after each operation on the vectors in the 


brackets, the components normal to V,, must be omitted. This follows from 
the definition of the relative operators. 
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and (V-N)? = (gr;- -‘N Ah — M2. 
Hence, choosing the remaining m—n—1 normals N,, to be orthogonal 
to N, and writing p»,,; = N-N,,;, we have, from (3.5), 
T+ M? = g''gikQ, 01,494 2 Poli olj> 


and hence, 


T = gg Q5;, Qh); Qn )+99 > boli Hols: 
Co 


From (3.21), we see that an alternative definition would be 
TM? — S(VN- NV), (3.23) 
where NV is the dyadic conjugate to VN, and S denotes the scalar. 
If m = n+1, we have 
| — gihgik ‘ 
Pe gg (Qj QQ, Qyy)- (3.3) 
The Gauss equation for V,, in V,,,, is 
— » ae Bay .49 
or 27, Q5,) ; Rijn 7 Rigs ~iy' yk Y’ iy we (3.31) 
where R, R’ are the curvature tensors of V,, V,,, respectively. 
Hence, (3.3) becomes 
T = eR—egiigiky By? 5y?), Rey 
e. by (2.3), 7’ = eR—e(a®—e€%€5)(ahy — oghe7), 
i.e. T = eR—eR'+2R,2E%¢8, (3.32) 
where Fg is the Ricci tensor for V,,, and R, R’ the scalar curvatures 


of V., V, 


n? m 


If V,,,, is an Einstein space, we have Rig = Ad,g, and (3.32) re- 


it 
duces to ‘ 
T = eR—e(n—1)d. (3.33) 
If the Ricci tensor vanishes in V,,,,, we have 
7 == €f. (3.34) 
In particular, this is satisfied if V,,,, is a flat space. For the case of 
a surface in V;, we have 
7 = § 2K, (3.35) 
where K is the total curvature. Weatherburn* showed that —2K 
was a generalization of the squared torsion, and the mean curvature 
a generalization of the curvature of a curve. 
We observe that, for a curve in hyper-space, 7' is the squared 
second curvature, and, from (3.33), this vanishes for a curve in Vj, 


the curvature tensor of a curve being zero. 


* Messenger of Math. 56 (1927), 173. 
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4. If the second fundamental tensors of V,, V,, in V, for the 

normals N,, are Q,,;;, Q,1.g, we have 
Oriap = —Ny0° Tp 
and Quigg = —Noe Ty = —Ny,e' rey: y? ;. 
Hence Qn = Qiapy iy? ;- (4.1) 
Also, from (1.41), we have 
Port = Noi Nuys = £5:°F a Nepy? ss 

i.e. Port = —Qriag Soi°Y? i. (4.11) 

If N, N’ are the principal normals and M, M’ the first curvatures 
of V,, in V,, and V, in V, ee we have 


MN = ¥ e,Qq)N 
and M'N’ = = 2 & Q,,N,, + 2% 0, N,. (4.2) 


Hence, the principal normal of V, in V,, is tangent to V,,, if 
o,.=0 (= eck44 ; ; (4.21) 
In this case, we see that 
N’ = N, M’ => M, Qi; — Q:;, (4.22 
and we may say that V,, is an asymptotic sub-space of V,,. These 
results are evidently true if V, is not necessarily a flat space. 
Writing p,,; = N-N, ., froma (3.22), (4.22) we see that, if 7’, T” are 
the second curvatures a V, in V,,, V,, in V,,, we have 
Pf <P Sansny (4.23) 
r 
i.e. from (4.11), 
T’—T = g'iy* ,y? gr > €- Qui Qrigss (4.24) 
id 
where N = £°r,. From (4.1), the equations (4.21) can be written 
giiy™ y? Quiag = 0. 
Hence, (4.24) can be written 
T’—T = gy; yP ; eves > e(Q, ay £2, p3—2, ap 2piy8)- (4.25) 
r 
This equation is true if V, is not a flat space. If V, is a flat space, 


the Gauss equations for Vi in V,, are 
> e,(Q, “ee 0, iaep 22, ya) saad Rrsyp- (4.26) 


Hence, from (2.2), (4.25) becomes 


(4. 
Tt” —-T = Rig EugB Rrsyp Eves 2 €5 Eqi*E ai? (4.27) 
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If m = n+1, we have €,.* = &*, and (4.27) reduces to 
T’—T = Rigéé. (4.28) 
Hence, from (3.32), we have 
T’ = eR—eR'+3Rg EE. (4.29) 
We observe that the direction of the vector é* is the principal 
normal direction of V, in V,, and 7” is the second curvature of V,, 
in V,. Also, R depends only on V,,. Hence 
Given a space V, in a flat space V,,, then all the spaces V,,,, which 
contain V,, and the principal normal direction of V,, in V,, at points of V,, 
have the same value for the expression 3Rig&¢?—eR’ at points of V,,, 
where Rig, R’ refer to V,,,,, and €* are the components in V,,,, of the 
principal normal of V,, in V,,. 
This could be written 
All spaces V,,,, of which V,, is an asymptotic sub-space have the 
same value for the above expression at points of V.,. 
For a curve on a surface in V3, we have n = 1, m = 2, p = 3, and 
Rig = —Ka,g where K is the total curvature. Hence we have 
: T’ = —K, 


R being zero for a curve. This is a well-known relation, for 7” is the 
squared torsion of the curve, and the curve is an asymptotic line on 
the surface. 

If V,, is any space, not necessarily flat, and if V,, is a sub-space 
ee 


geodesics in V,,. The necessary and sufficient conditions for this are 


is said to be totally geodesic if all geodesics in V,, are 


easily shown to be 


O,iap = 0 


r = 1,2,...,.p—m or = m—n-+1 
nf == 3,%....,m. 


From (4.26), we see that if V,, is a flat space, then all totally geodesic 


sub-spaces are flat spaces. 
If V,, is a totally geodesic sub-space of V,,, then (4.24) becomes 


T’—T = 0, (4.31) 


i.e. if V,, 1s a geodesic sub-space of a space V,, and if V,, is immersed 


»» then the second curvature of V,, in V,, is the second curvature 


in V, . 
of V,, in V,,. 


For the case n = 1, m = 2, p = 3, we have 7’ = 0 and hence 
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T’ = 0. If V, is Euclidean, V, is a plane, and the above theorem 
reduces to the property that the torsion of a plane curve is zero. 


5. Consider a system of o”-” spaces V, in V,,, such that one 
passes through each point of V,,. Then the normals N,, may be 


considered as functions of position in V,,. We have r; = r,y%,;, 
N,.; = Nogy8;, and hence, for the normal N,, = £,;°r,, we have 


Q,, = —9"t;- Noy = —giiy® ,y? ;r.°Nog, (5.1) 


7 omnes 


oj 


i.e. substituting for N,,, 

wit Q,, = —gtly* cy? 5 Eoia,p- (5.11) 
e have Ei, PE p= —é,i%p EwPEoia = —, ag ost 

where 9,)* = €y\*pEy 8 is the principal or geodesic curvature vector 


of the system of curves €,,* in V,,.. Hence, from (2.2), we have 


B 
O65 = — boa — 2 eu Map bo Me? (5.12) 
We have thus found the mean curvatures for the normals as functions 
of position of V,,. Obtaining the m—n mean curvatures in this way, 
we can now find the principal normal and the first curvature from 
the relations 
ME* = ¥ ¢, OES, M? = e> e,Q,?. (5.13) 
Having found the principal normal and first curvature, we can find 
the second curvature from the equation 


T+M? = g'[Ni-Niln = Fy iy AN ao Neplm (5.2) 
i.e., substituting from N = é°r,, we have 
T = —M?+(a* — Se, €,2¢,,P)a,5 £7, £9. (5.21) 
# 


If m= n+l, there is one normal £%, this being the principal 
normal, and we have 


‘---, um, (5.3) 
onl T = —M?+ (ax—ege¢P)a,5 £7, 2, (5.31) 
i.e. T = —M?+-aPa,; 7 , & —ee'x?, (5.32) 


where « is the principal curvature of the curves of congruence &%, 
and e’x? = ays nn, 1* = E% EB, 

From the above equations, we see that, given m—n directions, 
which we may take to be mutually orthogonal, at each point of V,,,, 
then, if these satisfy the conditions of integrability of the system 


of equations a,9£,,~dy? = 0 (o = 1,2,...,m—n), we can find the 
1 ap So) FY 
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principal normal and the curvatures of the spaces V,, defined as being 


normal to these m—n directions. 

Consider the case m = 3, n = 2, V, being referred to Cartesian 
axes. Let the given direction be (X, Y,Z), the components being 
functions of x, y, z. 

The condition of integrability of ¥ X dx = 0 is 

X(¥,—Z,)+ ¥(Z,—X;)+Z2(X,—Y,) = 0 (5.4) 
where X, eX /éx, ete. If this is satisfied, the above directions are. 
normal to a family of surfaces, and, from (5.3), (5.31), we see that 
the mean curvature K,,, and the total curvature, K, of these surfaces 


» 


are given by K.. = —(X,+Y,+%), (5.41) 


m 
2K = (X,+Y,+Z,)? + > [(XX,+ YX,4+ ZX;)?—(X?+4 X3+ X3)|. 
(5.42) 
As an example, consider 
= H 6 9 2) 9 
X > =- , re = L°--Yy* +2. 
> 
In this case, (5.41) reduces to K,, = —2/r, and (5.42) reduces to 
K = I1/r®. These are the results we should expect, for the. given 
directions are evidently normal to a family of concentric spheres. 








THE THEORY OF THE OSCILLATIONS OF CERTAIN 7 
STELLAR MODELS ABOUT THEIR STEADY- 
STATE CONFIGURATIONS (I) 


By L. E. LEFEVRE (Pasadena) 
[Received 7 March 1932] 


1. Introduction. The small oscillations about an equilibrium con- 
figuration of stellar models constructed according to certain hypo- 
theses as to physical structure, generation of internal energy, etc., 
have been considered by Eddington,* Jeans,+ L. H. Thomas,t and 
others. The investigations of Eddington apply a theory of ‘adiabatic’ 
pulsations of a particular model to the problem of the Cepheid 
variables. The oscillations of the interior regions of the star are sup- 
posed to proceed according to the ‘adiabatic’ law, pap’, and, sub- 
sequently, an estimation is made of the error involved in carrying 
.this assumption into the outer layers, by calculating the heat ‘leak- 
age’ there. A more satisfactory procedure is to make use of an 
equation of energy-exchanges which accurately expresses the balance 
at any point between the net gain of energy per unit mass due to 
generation of internal energy, less loss by radiation, and the various 
transformations of the different kinds of material energy by physical 
processes. Such an equation was used by Jeans,§ who obtained a 
third-order differential equation with respect to the time for the 
oscillations of certain models, and derived therefrom certain stability 
criteria. 

The outstanding difficulty of all existing theories of Cepheid varia- 
tion is their failure to account for the observed relation between the 
phases of the variations in displacement and luminosity. According 
to Eddington’s theory, the star throughout should be hottest when 
most compressed, whereas observational evidence shows that at the 


* M.N. 79 (1918), 2, 177. 

+ M.N. 85 (1925), 914; 87 (1927), 400, 720. 

+t M.N. 91 (1930), 122, and 619 (1931). These papers have appeared since 
the present investigation was begun. The fundamental equations of motion 
used by Thomas are similar to, but not identical with, those used here (see 
footnote, p. 305). The development of the argument differs widely from the 
present one, however, being directed toward a discussion of the nature of the 
stability criteria for stellar configurations. : 
§ e.g. Astronomy and Cosmogony (1928), Chap. IV. 
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surface the temperature is greatest when the velocity of approach 
of the surface towards the observer is greatest; that is to say, the 
luminosity there is in phase with dr/dt instead of r. 

In the general theory of small oscillations of systems possessing 
an infinite number of degrees of freedom, such as the vibrations of 
taut strings, organ-pipes, etc., the solution is obtained in terms of 
normal modes of oscillation, the most general motion of the system 
being a linear combination of such normal modes, whilst the funda- 
mental periods of the system are the roots of a certain equation 
which has to be determined from the boundary conditions of the 
problem. The boundary conditions in these problems are always such 
that they determine just one such equation, the ‘period equation’ 
of the system. Similarly, in atomic theory, the Wave Mechanics 
requires a study of the ‘eigenwerte’ and ‘eigenfunktionen’ of the 
differential equations. A rational theory of stellar oscillations, 
capable of application to variable stars, should aim at a differential 
equation, or system of differential equations, whose solution can be 
expressed in terms of normal modes of oscillation. No existing theory 
has established the existence of normal modes of oscillation of the 
radial pulsations of a spherical mass of gas evolving energy, nor has 
any placed the emphasis upon the influence of boundary conditions 
in determining the period or periods of such a system that is sug- 
gested by the analogy of the examples quoted. 

It is the object of the present paper to prove that, in accordance 
with general dynamical theory, there exists a definite “period equa- 
tion’ for certain particular stellar models. This is proved for con- 
figurations built on the ‘standard model’, both in the case of the 
completely gaseous (Eddington) model and also for certain ‘cored’ 
(Milne) models.* The existence of this equation is derived from a 
complete consideration of the ‘boundary’ conditions in each case. 
As, however, the differential equations of the motion are not soluble 


explicitly, the period equation is not given formal mathematical 


‘ 


expression. It further appears from these considerations that a 
definite relationship exists throughout the model between the phases 
of the variations in displacement and temperature. The determina- 
tion of this relationship would require the very troublesome numerical 
quadrature of the differential equations of the motion for certain 


* These models are chosen by way of illustration of the method only, and 
need not necessarily represent stars actually occurring in Nature. 
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suitable values of the disposable constants of the model. This has 
not yet been undertaken. 

The method adopted is similar to that of Jeans,* but the treatment 
differs from his in that we replace his single equation by a system 
of equations of lower order. 

In § 2 we briefly recapitulate for purposes of reference the equations 
of equilibrium of the gaseous layers of the Standard Model configura- 
tion.t The equations of motion for small oscillations of these gaseous 
layers are obtained in § 3. In § 4 we derive the system of equations 
for the most general possible harmonic oscillation of the gaseous 
layers as a whole. This system of equations (35)—(38) is the funda- 
mental one upon which the subsequent arguments depend. A very 
important feature of this system is the difference of sign between the 
last sets of terms on the right-hand sides of the last two of the equa- 
tions, suggesting the possibility of an accumulation of a difference 
in phase between the variations in displacement and temperature 
as we proceed through the model. It is pointed out in § 5, however, 
that these terms are only appreciable, in the case of stars of the 
masses, radii, temperatures, and periods of the Cepheids, in the layers 
lying close to the surface and so such a difference of phase is only 


appreciable in these layers. This is in agreement with Eddington’s 
conclusions.t The value of this change of phase in these layers 
depends upon the values of certain constants of integration which 
have to be appropriately determined from the boundary conditions, 


and owing to the complexity of the system of equations can only be 
determined by trial through the very troublesome numerical integra- 
tion of the system. It is significant to notice that when the extreme 
outside of the star is reached, the terms which were negligible in the 
interior become themselves the dominant terms in the equations (37) 
and (38), and so it is by no means certain that the phase change in 
the outer layers will be small, as Eddington, in the paper just quoted, 
has argued will be the case. 

In § 6 we consider the form of the solution of (35), (36), (37), (38) 
very close to the boundary, and determine the number of arbitrary 
constants left in the solution when the physical conditions of con- 
tinuity are inserted at the boundary. It is important to notice that 

* loc. cit. 

+ Eddington, M.N. 77 (1916), 16, and later papers ; Milne, M.N. 91 (1930), 4. 

{ M.N. 87 (1927), 539. See also J. J. M. Reesinck in M.N. 87 (1927), 414. 
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the boundary condition used by Eddington in his theory, namely, 
that 5P = 0 at the surface, is automatically satisfied, since we show 
that this is a necessary consequence of the circumstance that the 
amplitude of the displacement at the boundary must be finite. 

In §7 we state the boundary conditions to be satisfied at the 
centre of the completely gaseous standard model configuration, this 
case providing the simplest illustration of the arguments adopted, 
although it has become of less significance in the consideration of 
actual stars occurring in Nature, since the advent of Milne’s theory 
of stellar structure. It is established that these central boundary 
conditions are just sufficient to provide a determination of all the 
arbitrary constants left in the solution, together with a single equa- 
tion which may be interpreted as determining the period, or periods, 
of the system in terms of the fundamental constants of the model. 
One of thé conditions is, plainly, that there shall be no displacement 
at the centre. This is found to provide two equations between the 
constants left in the problem. It was suggested by Milne that a 
further condition to be satisfied at the centre is that there the star 
shall be hottest when most compressed, that is to say, the variations 
in displacement and temperature must. be in phase at the centre. 
This was assumed tacitly in Eddington’s theory and has not been 


explicitly stated previously as a condition of the required solution. 
This condition is fundamental in the present argument and provides the 
extra relation necessary for the determination of the period equation. 

In a subsequent paper we hope to consider extensions of the 


method to more complicated examples. 


2. The steady-state configuration. The equations of equili- 
brium of the standard model in which the opacity k and the rate 
of generation of energy per unit mass by internal processes € are 
constant, are already well known but are restated here for reference. 

Let M, L be the mass and luminosity, or total source strength, of 
the configuration respectively. Then if r denotes distance from the 
centre, p the gas pressure, p’ the radiation pressure, p the density, 
T the temperature, M(r) the mass within a sphere of radius r, L(r) 
the total net flow of radiant energy across such a sphere per second,* 

* The symbol L(r) is used in two senses. It can be taken to mean the total 
flow of energy across the sphere per second, as defined here, or the total amount 
of energy being generated within the sphere per second. In a steady state 
these are, of course, equal, but we shall be dealing with non-steady states in 
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G the gravitational constant, and c the velocity of light, we have for 


a steady state ep , ep’ _ _GM(r) | a) 
or ' or Alla 
ep’ ——ikL(r) (2) 


or 4ncr? Ps 


| ly) _._ 
and Mr) oe ae (3) 


From these (assuming as boundary conditions p = 0, p’ = 0) is de- 
rived the familiar relation 


p= 


kL 
} inate = ]J— 
bg p dncGM’ 
a constant of the model. 
In the outer (gaseous) regions of the star the equation of state 


p= Os 
pe 


holds, and from this and the relation 
dM(r) 
ed 
dr ide 
is derived the Emden ‘polytropic’ equation 
1 d/,,d0 
— —|@—] = —#@ 6 
zig a) - —* rm 
in which the variables have been changed according to the relations 
7 = A, (7) 


_ '(4e-F 3 
r= ala) Bp Geaay? Bi 


and es = 
The ‘discriminant’ of the model is 
a 1—f/R\* a I 
4 pt in draG?M’ 


a known function of L, M, and k. 





which this is not necessarily true, and so it is necessary to keep the distinction 
clear and use the symbol only in the sense defined here. 
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If A is so chosen that €,, the value of € corresponding to r = r, 


the radius of the star, is equal to unity, then 


= R\ (1 -B)} ] 1 : 
r= (7) BI (10) 


2. 
bra)! r, 


B 


and the boundary value of d6/dé is given by 


dé ] 
(2), ree “- 


If L has the critical value L, of Milne’s theory the configuration 
is wholly gaseous and of the type investigated in Eddington’s theory. 
Otherwise the configuration is either of the ‘collapsed’ type (ZL < Lp) 
or the ‘centrally-condensed’ type (L > L,) of Milne’s theory. We 
shall take note of this distinction later when considering the boundary 
conditions for the different models to be discussed. For the present 
we confine our attention to the derivation of the equations of motion 
of the gaseous layers of the star. 


3. The equations of small motion of the gaseous layers. We 
fix attention on a particular element of material which in the equili- 
brium configuration is distant r, from the centre. The physical 
conditions of this element in the equilibrium configuration will be 
represented by symbols with zero suffixes attached. In the small 
motion to be considered we put 

r= 1)+6r 

T = T,+8T 

p = potdp 
etc., where 57, 57’, 5p, etc., are small compared with 79, 7), po, etc., 
and are functions of 7, and of the time ¢. 

The equations to be satisfied throughout the small motion are: 

(a) a hydrodynamical equation of motion 

Du __ = 1 &(p+p’) GMi(r) 


Dt p or r? 


where w is the outward radial velocity and D/Dt denotes differentia- 


(13) 


> 


tion following the motion; 

(b) an equation of energy-exchanges, expressing the fact that the 
rate of gain of energy per unit volume of any portion of the material, 
less the rate of loss of energy by radiation, is equal to the rate of 
gain of energy measured in terms of the physical properties of the 
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material: such an equation has been derived in a general form by 
Milne:* the form appropriate to radial motion of a spherical con- 
figuration is 


; 2 DE 
cm eg NR ae ee 
€p cfr ()}= + 
where F(r) is the net outward flux per unit area across a sphere of 
radius r, E the density of radiant energy per unit volume, and c, the 
specific heat per unit volume of the material; 

(c) an equation of continuity. This is obtained by expressing the 
fact that the mass within a small shell of radii r, and r,+-Arg remains 
unaltered in the disturbed configuration, and is easily found to be 

dp , Wr adr 
7 ee = (15) 
Po % 

In substituting the expressions (12) into equations (13) and (14), 
the following lemma, which is capable of elementary proof, is of 
considerable use: 


Lemma. If ¢ is any function of position in the sphere and ¢ = dy 
when r = 1%, in the equilibrium configuration, and ¢ = $)+8¢ when 
r = ro+6r, in the disturbed configuration, then the space gradient of 
d, 0b/er, in the disturbed configuration is equal to 


(a) ) 


neglecting quantities of the second and higher orders of smallness. 


We can therefore write (13), to the first order of small quantities, 





_ i (p+P')\ ()_ 2r) + Ses P)) _SO—*] 
Po Po or 0 org org r Yo ‘ 


since M(r) = [M(r)],) always. 
Using (1) and (15) this becomes 


d*(5r) _ — > (Ae?) br 1 (p+p') (16) 
Po 0 


dt erg T Po 

* Quart. J. of Math. (Oxford), 1 (1930), 1. Other forms of this equation 
have been given by H. Vogt in A.N. 5545 (1928), 232, and L. H. Thomas in 
Quart. J. of Math. (Oxford), 1 (1930), 239. 


3695.3 x 
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Similarly, to the first order of small quantities, (14) can be written 


Sr\ COL (r a(Sr dds L(r)) 
(€9+ 5e)(py+Sp) ie zal! = = )K( 2 ) ( " = ) ” = | ie 
, rr ¥ 0 0 “- 


ah Yo 
/ 


of d(3p = (F) dsr) | Po Ty U(5p') (p+p’) d(dp) 


dt or Jo dt = = 4p dt Po ae 
since H = 3p’, p’  T4, and F(r) = L(r)/4ar?. 

The coefficient of (@L(r)/ér), on the left-hand side of (17) is, to the 
first order of small quantities, 

"ys 25r f. or) 
; ro cro 

and this, by (15), is equal to , 
tao 


Po 


eL(r) ie ee) 
or Jo of or jo 


= €9 47717 po- 
Hence the left-hand side of (17) reduces, to the first order, to 
*) 1 ed L(r) 


4nr, Oro 


> 


€o Po7 <8 pude—coPa T 


Po 
1 aL(r) 
4nr, Oro : 


DT) 2 
Further, Cotte To: Po_. 
y—1 


that is, to py de — 


where y= c,/¢,, 
the specific heat ratio for the matter alone. Hence 
Cy po Ly i B 
4, ~ M(y—1)(1—B) 
from (4). 
Hence we can now write (17) in the form 


12(y—1)(1—B)-+B dp’) 

4(y—1)(1—8) dt 
— 3(° d@r) | (pt+p )o dp) . 56 L OL (7) (18) 

or }, dt Po dt 4nr= Org 

Owing to our lack of knowledge of the nature of the processes by 
which energy is liberated sub-atomically, we have no information as 
to the manner in which S¢ will vary during the motion. We must 
here introduce an arbitrary hypothesis as to this variation. The 
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simplest suitable hypothesis to make is that 5« = 0. The introduction 
of this hypothesis is equivalent to an averaging process for the energy 
liberation throughout the duration of the motion, necessitated by 
our lack of knowledge of the actual variation at each instant. It 
should be clearly understood that the hypothesis is quite distinct 
from that made in taking « to be constant throughout the equilibrium 
configuration, which is equivalent to an averaging of the energy 
generation throughout the star. Each hypothesis represents an 
averaging process, the one in time, the other in space. 

Our equations are then, finally, 


@Gr)__ 4 0(p+p')8r_1A8(P+P') aw 


dt? p oO rp or 
12(y—1)(1—B)+B dp’) _ 3ap' dr) _ p+p' dbp) 1 e@Li(r) 


4(y—1)(1—B) dt  o& dt p dt 4ar® or ’ 
(20) 
5 2r , a(6r 

= — (74, (21) 
p r or 

in which we have dropped all the suffixes in the coefficients of the 
variational terms, it being understood henceforward that all such 
coefficients refer to the equilibrium values of the physical quantities 


concerned. 


4. Conditions for harmonic oscillations of the gaseous 
layers. We seek the most general harmonic oscillation of the 
gaseous layers in which each element performs harmonic oscillations 
of period 27/n. Such an oscillation will be represented by the 
expressions 


dr = A(r)cos nt + B(r)sin nt ) 
dp’ = C(r)cosnt + D(r)sinnt |’ 
where the functions A(r), B(r) are small compared with r, and C(r), 
D(r) are small compared with p’. (A knowledge of these functions 
for the gaseous layers would determine the relative phases of the 
variations in displacement and net flux at any point; in particular, 
at the surface.) 
Then, by (21), 


ee (74 +22 cose + 
p 


(22) 


r dr 
oT 1 dp’ 


, £2 


(2497) sin nt}, (23) 


x2 
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sT ce } 
— = -——{Ccosnt + Dsinnt}. 
T 4p 
dp. . _ 8p 67 
Pp p ve 
— C cos nt + Dsin nt 

4(1—B) p 
using (4). 

Hence 


Hence 


—3B 
rh 


d(p+p’) = C cos nt + Dsin nt)— 


. qT 9 T) 
— AP (= Pi =) cos nt += 2B +5) sin nt}. (25) 


Also, from (2), 


dL(r) - 26r _#.. op’ i; p’ 
Lir) oT p or }/ er 
a ce ye Es 


7 7 7 r z ? 
r or or or or or 


using (21) and the lemma of § 3. Hence 


we) (C4 dC Peo cnt + (43 5 0B [2 in nt (26) 


Lir) r ‘dr r  adr| oo 

Substituting these expressions into equations (19) and (20), and 
equating the coefficients of cosnt and sinnt separately to zero, we 
find that the functions A(r), B(r), C(r), Dir) must satisfy the system 


of differential equations: 


2d — 1 &' 7, 1 4-38 an 
TB pr or * 4(1—B) dr 1—Bér 
7 4 1ép',,1f 4-38 dB pB é 
B+- = a 
sia ~ 18 pr or “lai dr 1—Béor 
12y—1)(1—P)+B 6 
4(y—1)(1—B) 
> op’ y np’ (2A dd c @ {r? 4B ap’ 
= A= —— ass ec eae s 
~S rte dr)” kr ér\p\ r ar 
1%X(y—1)(1—B) +B 5 
4(y—1)(1—B) 
Op’ ~ np" (2 a + c é (a(t op’ dd 


B : 
r kr? ér | p 


as By 
° Or 1—f 


r Or dr 
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We express these equations in terms of the Emden variables 
defined in § 2. We write 


_1/R\(1—p)F 1 
A= 3) a aaa 4 





or, by (10), A(r) = r,A(é), (31) 
where r, is the radius of the star, and, similarly, 

Bir) = r,Bé). (31’) 
Also, put O(r) = sadr4C(€), (32) 
and Dir) = sar*D(€). (32’) 


We further introduce new symbols w? and A defined by 
ee: a 


~ 4 (Rjp) 1-8" 
or, in terms of the discriminant of the equilibrium model, 


WwW 





me! = - . 
~~ 1-8 eur’ _ 
e 1 1—B (R/p) 
é 1 / SS eek cen ee se : 
unc A in Bp jae 
1 ly 
pene 34 
i “= Ramee (34) 


(We observe that 2w? and A/n are dimensionless numbers.) 
Then, using the relations 
p’ = far*4 
1 a Bp 7 
and , ome (Riu) 1B 
the equations (27), (28), (29), (30) become (dashes denoting dif- 
ferentiations with respect to &) 


so — (oes) 


4(1—) BE z 

Se ois oe atin seh) 
12(y—1)(1—B) +B o 

4(y—1)(1—B) 


. 64 2A \ Alé@ 47°B. D’ 
= 3074 ae t4)-naagl(g tap 
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12(y—1)(1—B) +B p 
4(y—1)(1—B) 
: 64 2B Al feo 40’A , C’\) 
= 366’ B_____—_ + B’)+-.=— — . & 
wile : +; r BE" eC; ro =) i 
The variations in the various physical quantities involved are 
found to be, in virtue of equations (22), (23), (25), (26), and the rela- 
tions between the new and old variables, 
dr = r,(A cosnt + Bsin nt), (39) 
-B GM2C 
4 re 
_ MCG (/24 
4a} | é 


op —(C cos nt + Dsin nt), (40) 


Sp = - +4’) cost + (= +2')sinad, (41) 


GM2C 


d(p+p’) = — ja 38)(C cos nt + Dsin nt)— 


— pot (= +A ) cos nt + (= + B') sin nt | (#3) 


sL(r) = — A -|(FF+ 48 fa) 008 nt + (“e+ =) in ntl. (43) 


l6zrt 


ie J 

The equations (35), (36), (37), (38) are four linear differential equa- 
tions for the functions A, B, C, D. The solution of the problem of 
the small oscillations of the gaseous layers of a standard model con- 
figuration consists in the determination of a solution of these equa- 
tions, which satisfies appropriate boundary conditions according to 
the model under consideration. These conditions for certain parti- 
cular models are considered in later paragraphs. Such a solution 
being determined, the physical changes of the configuration during 
the motion are completely determined by equations (39), (40), (41), 


(42), (43). 


5. Simplification of the equations in the interior of the star. 
The last sets of terms on the right-hand sides of (37) and (38) have, 
it should be noticed, opposite signs. These terms arise from the 
term involving 5Z(r) in (30), and hence will be for convenience 
referred to as the ‘flow’ terms. The difference of sign in the flow 
terms in equations (37) and (38) suggests the possibility of an accu- 
mulation of phase difference between the variations in displacem ent 
and luminosity as we proceed through the star. But it must be 


4a | é 
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* noticed that for a star whose internal temperature is of the order of 
several million degrees, oscillating with a period of the order of a few 
days, the quantity A/n, defined by (34), is very small. So long as 
6° is large in comparison with A/n, which is true for a very great 
part of the star, and fails only for a very small fringe at the surface, 
the flow terms in (37) and (38) are negligible. When we neglect 
these terms, the system of equations falls naturally into two pairs, 
one involving only the functions A and C, and the other involving 
only B and D, and the latter pair is the same as the former except 
for the replacement of A by B and of C by D. In these circumstances 
A/C = B/D for all values of £; that is to say, there is no change of 
phase between Sr and dp’ in these regions of the star. Any such 
change is due to the effect of the flow terms in (37) and (38) and 
will therefore occur only in the outer fringe of the star when these 
flow terms become significant. This is in agreement with the results 
of Eddington’s theory of ‘adiabatic’ pulsations. 


6. Form of the solution of the equations near the surface. 
The system of equations (35), (36), (37), (38) consists of four linear 
homogeneous equations, each of the second order. The most general 
solution will therefore contain eight arbitrary constants of integra- 
tion. Eight ‘boundary’ conditions will be required to determine 
these, and one further will be necessary to provide a period equation 
for the system; this latter equation is, then, the eliminant of the 
eight constants from nine equations. It is not at once obvious, how- 
ever, how many conditions are represented by the circumstance that 
the amplitude of the oscillation must remain finite at the boundary, 
i.e. that the functions A and B in (35), (36), (37), (38) must be con- 
tinuous at the boundary. It is therefore first necessary to investigate 
the most general continuous form of the solution of these equations 
very close to the boundary, ¢ = 1, with a view to discovering the 
degree of indeterminacy left in the solution when the condition of 
continuity of the amplitude at the surface has been satisfied. 

Let 7 = 1—&, so that 7 is small near the surface. 

Then it can easily be shown that the form of 6, the solution of 


(6), near the surface, is 


6 = = {n+ +9°+ 4°+0(n°)}. 


Then, if A and B are continuous at the surface, the right-hand 
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sides of equations (35) and (36) are at least O(n*) very close to the 
surface. 

Hence C’ and D’ are at least O(7*) very close to the surface; hence 
C and D are at least O(n*) very close to the surface. 

(It follows from (42) that, at the boundary, 5P is of necessity zero, 
if A and B are continuous there. The condition 5P = 0 cannot 
therefore be used, as Eddington uses it, as a boundary condition to 
determine the appropriate value of n? for a physically suitable solu- 
tion, since it is a logical consequence of the equations of motion 
themselves, together with the physical property of continuity of the 
solution at the boundary.) 

The left-hand sides of (37) and (38) and the right-hand sides with 
the exception of the flow terms are all therefore O(7*) at least. Hence, 
sufficiently close to the boundary, these flow terms become dominant, 
however small A/n may be. (A/n = 0 implies n infinite.) 

Hence, sufficiently close to the boundary, we can neglect all but 
the flow terms in (37) and (38), so that 

, Vet 3) + = 0, 
{ é 6 }) 


and . 46" A, = 
‘ 6 + a ie 


! ) 

Hence 4. = 
40'A OC’ 
& ria © 


where k, and k, are arbitrary constants. 


and 


Further, the dominant terms in (35), (36), namely, those in 6°, give 


4-38 C’_ 4 O\ 4, B 924, 4 


4—38 D’ ee 2 ae 
é —__ __ = | w*n* — £.. é 
und 70) Hf * a n , ad +B) 


Hence, sufficiently close to the boundary of the star, the function 
A must satisfy the equation 


ama (rpg) tee" (e +4) acne 
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and the same is true of B, k, replacing k,. Hence, 
1—B 0’ 4— 38 k 
'A'+ ("FP nt +s “A= Ka 44 
p g)°— Be - 
+p , (1—B 0’ 4— 38k 
and 0B +( : ntot—2) B= : = 
B g 4p 


In the region under consideration the form of @, given by the 
solution of (6), is 


(45) 


* 0 = (nb att atta tly 1°), 


where € = 1—». 
so that 7 is very small indeed. Thus, to a close approximation, 
ocak, 
Ct 1-7’ 
Li-g 
= s 
ee 1 
ce 
Hence we may rewrite (44) and (45) in the form 
4 =F 
48 
a ee ~, 4—38 . 
B’—|— ntoCige + 2) B = —C} kp. (45’ 
.. : rd 
Integrating these equations we find that the most general con- 


tinuous form of the functions A and B very close to the surface of 
the star can be written 


: = 6% = | get atl, 


= 


Then 


(44’) 


A (PE ntatOige )4 = = —C- 


= £-le- {KC at “ p | bere ae|, 
where K, and Kp, are undetermined constants and 


1—p 2 27% 
= —" n®w*Ct., 
B 


The solution thus contains four arbitrary constants, ky, kp, K4, 
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K,. We have as yet made no stipulation concerning our origin of 
time. If we so choose this origin as to fix the variation in energy- 
output at the surface in value and phase, we at once assign values 
to k, and ky. Thus, if we choose the origin of time so that the 
variation of L(r) at the surface, i.e. the observed variation in bright- 
ness, is a pure cosnt variation, we can put k, = 0 and assign a 
definite value to kp, according to the amplitude of the variation in 
brightness. Suppose, then, we choose k, = 0 and kp = 1. Then the 


solution contains two remaining undetermined constants, K, and K p. 
Knowledge of the ratio of these two constants would give us the 
phase of the variation in displacement at the boundary. For when 


San Rew, 
and B= e°K,, 

since we may, without loss of generality, assume the additive con- 
stants for the indefinite integrals to be such as to make these integrals 
vanish for € = 1. But the ratio cannot be determined by considera- 
tion of the surface layers alone. The constants K, and Kp, can only 
be evaluated by finding which solution of (44) and (45) can be fitted 
on to the solutions of the equation in the interior regions where the 
present approximations break down. 

We have now considered the solutions of the general system of 
differential equations for the motion of the sphere in two different 
regions—in the interior regions where the flow terms are negligible, 
and in the surface layers where the flow terms are dominant. The 
former regions extend throughout almost the whole sphere; the latter 
occupy but a small fringe at the boundary. Between these two 
regions there lies a transition region where the flow terms are com- 
parable in magnitude with the other terms and where the equations 
of motion do not admit of simplification. The intractable nature of 
the equations in this region leaves a gap which cannot be bridged 
mathematically and makes it difficult to link up the appropriate 
solution for the interior with that for the surface. We are thus unable 
to determine quantitatively the phase relation between the variations 
in displacement and temperature at the surface. As we have already 
remarked, these variations are in phase in the interior or ‘adiabatic’ 
region. It is in the transitional region that a phase difference begins 
to accumulate, but by how much we cannot say. K, and Kp, can 
only be determined by continuing the solution for the interior 
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through the transitional layer to the surface region, and though this 
is possible theoretically, the nature of the equations makes it im- 
practicable. 

The region in which the approximations (44) and (45) are legiti- 
mate depends upon A/n. The necessary criterion is that @ shall be 
negligible in comparison with A/n. Now for spheres of the masses,_ 
radii, and periods of the Cepheids, the value of A/n is about 10-*. 
Hence @ must be negligible in comparison with 10-7. Thus, for 
example, to four significant figures the approximations are legitimate 
only if 6 < 10~*, and this is true for so minute a fringe at the edge 
of the star that the approximations described above are of no value 
in the actual solution of the equations, although, as we have seen, 
they indicate the degree of indeterminacy left in the problem when 
the surface conditions have been chosen appropriately. 


7. Internal boundary conditions for a gaseous sphere. As 
already indicated, the determination of K, and Ky requires the 
choice of an appropriate solution in the interior region. It therefore 
remains to consider what internal boundary conditions must be 
satisfied for certain models. In the present paper we confine our 
attention to the simplest model available—the completely gaseous 


sphere (Eddington model).* In Part II of this paper more com- 
plicated models, based on the recent researches of Milne, will be 


considered. 

In the present case the conditions to be considered are those at 
the centre, € = 0. 

In the first place it is obvious that two conditions to be fulfilled 


there are A(0) = 0, (46) 
and B(0) = 9, (47) 


since the centre must remain stationary. 

A further condition that must be satisfied at the centre is that 
the phases of the variations in displacement and temperature as the 
centre is approached will be such that the instant of greatest com- 
pression there will coincide with that of highest temperature.t There 
is, of course, no direct evidence for this, but physical intuition is 


* Jeans has claimed in the series of papers quoted in the Introduction to 
this paper that this model is vibrationally unstable. 
+ Suggested to me by Professor Milne. 
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strongly in favour of its validity in the case of matter at rest, 
generating energy at a constant rate (since we have assumed that 
de = 0), while subject to periodic variations of compression. 
Now, dr oc A cosnt + Bsin nt, 
and dT x Ceosnt + Dsin nt. 
d(sr) 
dt 
d(5T) 
dt 
Hence, greatest and least compressions occur at the centre when 


tan nt = B/A. 


Hence oc n(—A sinnt + Beosnt), 


and oc n(—C sin nt + Decosnt). 


Similarly, greatest and least temperatures occur when 
tan nt D/C. 


Hence the above condition requires that 


(48) 


The equations (46), (47), (48) constitute a set of three relations of 

the type = Fin, 8, y,01,K, Kz) =0 (i =1,2,3). 
From them K ,, K, can be eliminated, leaving a relation of the type 
d(n, B, y, C!) = 0. (49) 


where ¢ represents a function not necessarily single-valued. 
This is the period equation whose existence we desired to establish, 


appropriate solutions for the system being obtained only for such 


values of n as satisfy this relation. If this equation has its roots in 
n all real, the only possible motions are stable ones, but if imaginary 
roots in » are possible, unstable motions may occur and the model 
is unstable vibrationally. As (49) cannot be written in explicit form 
from the present analysis, no examination of the nature of its roots 
is here possible. 

Further, when an appropriate value of x has been chosen according 
to (49), any two of the equations (46), (47), (48) determine corre- 
sponding values of A, and K,. This determines the phase relation 
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between the variations in displacement and temperature (and there- 
fore also in luminosity, from (43)) at the boundary. It is therefore 
established that a definite phase relation does exist at the boundary 
between these variations. As stated earlier, this could only be deter- 
mined through the long and excessively complicated numerical 
quadrature of the fundamental system of equations for specimen 
values of the constants involved (8 and y). 

Note. In the course of preparation of the present paper for publica- 
tion there has appeared a thorough investigation by 8. Rosseland* of 
the nature of the stability criteria for gaseous spheres subject to 
radial pulsations. The procedure employed is essentially the same as 
that adopted here. (The equation of energy-change used by Rosse- 
land differs from Milne’s equation, being identical with that used by 
both Jeans and Vogt in their earlier investigations, quoted at the 
commencement of this paper.) While Rosseland’s results are, in the 
main, directed towards stability considerations, and, for stars built 
on the standard model, corroborate the conclusions of Jeans, obtained 
by less elegant methods, his investigations of the atmospheric pro- 
blem lead to results essentially the same as those of the present paper 
as regards the extent of the ‘adiabatic’ region of the star and the 
accumulation of phase difference between the variations in displace- 
ment and temperature in the transitional region just below the sur- 
face when these equations of motion do not admit of simplification. 

The development of Rosseland’s paper follows the :ines of the 
general dynamical theory of oscillating systems and goes far towards 
supplying the need, mentioned in the Introduction to this paper, of 
a rational theory of stellar oscillations. Thus, for example, the ortho- 
gonality property of the oscillational functions for the amplitudes is 
derived for certain systems. Also an estimate is made of the order 
of magnitude of the period of the fundamental mode of oscillation 
of the star. The establishment of a period relation by the detailed 
consideration of boundary conditions, which is the main feature of 
the present paper, does not, however, form a part of Rosseland’s 
investigation. 


* ‘On the stability of gaseous stars’: Univ. Obs. Oslo Pub. No. 1 (1931). 














’ A’ RELATION BETWEEN HYPERGEOMETRIC 


SERIES 
By J. L. BURCHNALL (Durham) 
[Received 9 November 1932] 


Mr. H. B. C. Darwtne has recently* given the analogues, for func- 
tions of higher order, of the formula 


Cn Bes oe oe oe 
. 9 4 


¢ 2 2—c c 


The machinery he employs increases in complexity as the order of 
the functions increases and does not readily establish the general 
result suggested by the particular instances. I venture therefore to 
give a general proof from which it will appear that the theory of 
adjoint operators indicates the existence of relations of the type in 
question and that a simple partial-fraction theorem suffices to deter- 
mine their precise form. 

In the interests of symmetry I slightly modify the conventional 
notation. Let c,, ¢.,..., ¢, be numbers incongruent to modulus unity 
of which one, say ¢,, is unity. Then I write 

r( @1, Ug...) Bn 2 )= ee = Glace. i) 
Se ere Cy, Cqyeeey Sms 

Let 5 denote ad/dx, and notice that the operator adjoint to +a 
is 6+1—a. 

Consider the equation 


TI (8+0,—ly = «TI (+a,)y. (2) 


r=1 r=1 


[ts adjoint is II (8+2—b,)w = [[ (6+1—a,)aw. 


r=1 r=1 


If z = aw, then z satisfies 
TI (8+2—6,—1)z = « TJ (8+1—a,)z. (3) 
r=1 r=1 

Linearly independent solutions of (2), (3) are, for r = 1, 2,...,n, 


y, = «1-+G(a+1—b,,b+1—b,,x) = ai y A,” 
D m=0 ; (4) 
z, = x-1G(b,—a,b,—b+1,2) = a ¥ Ba 


m=0 


* Proc. London Math. Soc. (2), 34 (1932), 323. 
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where 
is =I [(a, +i b ml (b + 1—b baad 
(5) 


B — II [(6,— s)m/ (6, —b,+ 1) mn] 
Now, by the general theory of the adjoint,* there exists an identity 
of the form > Ae Yr% = 0, 


and consideration of the various terms in the neighbourhood of the 
origin makes it clear that the identity is of the form 


Yb Ye = 0. (6) 
More precisely I show that, if f(t) = II (t—b,), then 
r=1 


2 [yr *if'(b,)] = 0. 
It is evident that (7) is established if, for every k, 


E[( 5 4rn-m Brn)/f’Go)] = 0. 


To obtain (8) write g(t) = Il (t—a,) 
r=1 


and consider the ‘key-function’ 
g(t—1)g(t—2)...g(t—k) 1“ 

FOft—1) ...f(t—*) 

The analysis has proceeded on the assumption that the b, are incon- 
gruent to modulus unity, so the denominator of H(t) has no repeated 
roots. Moreover, the degree of the denominator exceeds that of the 
numerator by at least two, since n is assumed greater than unity. 
Thus the sum of the residues of H(t) at its various poles is zero. 

Now the residue at t = b,+-m (m < k) is 


m—1 


TI[II (,+1— ad| TTT (b,-I— a] 
{i [I] +! of’ TE [IT 1-00] 


H(t) = 





‘ )n( a,+1— Dy) i oY 
gs I lo, = bt Vlg F1—Byy “lroy 
which is Uacaill m/f '(b,), and, on summing over all the poles, 


(8) is established. 


* See, e.g., Darboux, Théorie générale des surfaces, ii, p. 103. 
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It is to be noted that the theory of the adjoint* requires also the 
following relations between the y, and the ith derivatives z\ of the 


namely, 
> [9-4 /f'(b,)] = 0 (¢<n—l) \ 
= Cx1-"(1—2)-1 (i = n—1) J 


where C is a numerical constant. On picking out the coefficient of 
x-" on the left of (10), we have 


nr 


» (10) 


, XO (b,—1)(6,—2)...(b,-—n +1) _ 
o=)> f'(b,) ica 


by a well-known theorem in partial fractions. 


(11) 


r=1 


The ‘key-function’ corresponding to (10) is 
H,(t) = (t—1)(t—2)...(t—i)H (0). 
Since the derivative of a hypergeometric series is a multiple of a 
hypergeometric series, the relations (10) lead to a great variety of 
relations between such series. The simplest is, in the usual notation, 


(1—a)(1—6) P (" ae eg i 
(1—c)(2—c) c led 


2—c c+1 


a(t 1—c,b+1—c; ") steal POP 4 


c(c—1) 
aed rae 5 a 


2—c c 
= «1(1—2)-1. 


* Darboux, loc. cit. 











HEFFER’S FOR ALL BOOKS 
Scientific Books and Journals of all kinds, pete and Second-hand, 
English and Foreign. Heffer’s are always repared to purchase 
complete sets or ~~ runs of the Quarterly 7. Ty Journal of M 
Catalogue 386, Mathematics, Physics, Chemistry, and 
Publications of Pace, Societies, free on request. 


W. HEFFER & SONS, LTD., CAMBRIDGE, ENGLAND | 





OXFORD BOOKS 


WAVE MECHANICS 
The Elementary Theory 
By J. FRENKEL 
20/- net 


THE THEORY OF ELECTRIC 
AND MAGNETIC SUSCEPTIBILITIES 
By J. H. VAN VLECK 
30/- net 
ai ia Cae eae cine looney amend UO denied wots a 
excellent reviews are given of experimental results. The book is a most ¢ valuable 


contribution to the weed of magnetic and electric susceptibilities, and to Quantum- 
Mechanical Theory. . . .”—Nature 


THE PRINCIPLES OF 
QUANTUM MECHANICS 
By P. A. M. DIRAC 


17/6 net 
“|. . There can be no doubt that his work ranks as one of the high achievements 


of contemporary physics. . . .”—Nature 
THE CORRESPONDENCE AND PAPERS 
OF EDMOND HALLEY 
By E. F. MacPIKE 


.21/- net 


. This beautifully printed volume gives an admirable Ar into the scientific 
activities of one of the most famous of English astronomers. —Cambridge Review 


OXFORD UNIVERSITY PRESS 



































DEIGHTON, BELL « CO., LTD. 
13 TRINITY STREET, CAMBRIDGE 


UNIVERSITY BOOKSELLERS 
ENGLISH AND FOREIGN 


We supply 
NEW AND SECOND-HAND TEXT-BOOKS in all Subjects 
FOREIGN BOOKS in all Languages 
The LATEST PUBLICATIONS in all Branches of Literature 

OLD AND RARE BOOKS, &c. 
We buy 
LIBRARIES Or SMALLER COLLECTIONS OF SCIENTIFIC BOOKS 
JOURNALS, PUBLICATIONS OF LEARNED SOCIETIES, &c. 
English and Foreign 





Orders by post receive prompt and expert attention 
Catalogues issued regularly 


Telephone: CAMBRIDGE 239 ESTABLISHED 1700 














a 




















Published in 1932 


FOUNDATIONS OF POINT SET THEORY. By 
R. L. Moore. American Mathematical Society 
Colloquium Publications. Volume XIII. 7 + 486 pp. 

$5.00 


DIFFERENTIAL EQUATIONS FROM THE 
ALGEBRAIC STANDPOINT. By J. F. Rirrt. 
American Mathematical Society Colloquium Publi- 
cations. Volume XIV. 10 + 172 pp. $2.50 


All volumes of the Colloquium Series may be ordered from the 

American Mathematical Society, 501 West 116th Street, New 

York, or from Bowes & Bowes, 1 Trinity Street, Cambridge, 
England. 





Printed in Great Britain at the University Press, Oxford, by John Johnson, Printer to the University 


ot 


\t 





























